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ABSTRACT: Humic Substance (HS) samples were collected from twenty-five sites along the 
Suwannee River. These sites included the Okeefenokee Swamp, several springs and tributaries, and the 
delta where the river empties into the Gulf of Mexico. Fourier Transform Infrared Spectroscopy (FT-IR) is 
used to categorize the relative abundance of the functional groups of these samples. In salt-water samples 
we identified a relatively intense C-Cl vibration indicating the chlorination of the humic substance. We 
also measured the vibrational spectra of humic acid with several concentrations of HNO3, KOH and 
NazP207 over twenty-eight days to optimize our extraction procedure by minimizing condensation 
reactions. 


Key Words: Humic, fulvic, humic substances, infrared, IR, FT-IR, Suwannee 
River 


Humic substances (HS), the product of plant and animal decay, are divided into 
three classifications: fulvic acid (FA), humic acid (HA), and humin. FA, the lowest 
molecular weight component, is characterized by a yellow to tan color in solution 
and is soluble in water at all pH’s. HA is insoluble below a pH of two and is brown 
to dark brown in color in the aqueous phase. Humin is insoluble in water at all 
values of pH, is black in color, and is found at the bottom of rivers, streams, oceans, 
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and lakes. The carbon content of FA is lower than HA and humin, and the mass 
percent of oxygen decreases with the progression from fulvic acid to humic acid and 
onto humin (Pettersson et al., 1994; Schulten, 1996). HS are macromolecules that 
aggregate and are comprised of cellulose, lignin, amino acids, peptides, sugars and 
contain significant quantities of the functional groups phenols, amines, carbonyls, 
and carboxylates. A humic acid model proposed by Schulten (1996) has the 
chemical formula of C39gH33sOo9Ns5 giving mass percentages of C 66.69%, O 
25.96%, 6.09% H, and N 1.26%. The International Humic Substances Society 
(IHSS) uses humic acid and fulvic acid from the Suwannee River at Fargo, GA as its 
standard (Averett et al., 1994; Manning et al., 2000; Ouatinane et al., 2000; Rasyid 
et al., 1992). The molar mass of FA from the Suwannee River (Fargo, GA) was 
reported to be 1000—1500 amu and molar mass of HA from Fargo was reported to be 
5000-10000 amu (Thurman et al., 1982). When Ouatmane and coworkers (2000) 
used FTIR to study HA, they found several reoccurring spectral features including 
an O-H group, an aliphatic CH stretch, the presence of C=C, C=O, and/or C=O of 
conjugated ketones, quinones and amides, spectral evidence indicated C-H 
deformation of CH, or CH groups, and correlated spectral features below 1290 
cm! with aromatic C=C, COOH and CHs3, C-O of aryl ether, and C-O and OH of 
COOH groups. Their list does not represent a complete list of potential vibrations 
that could be found in Suwannee River Humic Substances (SRHS), such CN or C-Cl 
vibrations so a more complete table of functional groups and frequencies that can be 
observed in HS are presented (Table 1). 

Dupuy and Douay (2001) studied the behavior of heavy metals in soils that 
complex humics by infrared spectroscopy. Klucakova and coworkers (2000) 
characterized Humic Acids (HA) and fulvic acids (FA) by their ability to form 
complexes with inorganic and organic molecules capable of causing the effect of 
cross-linking by FT-IR. Lynch and Smith (1992) reviewed inconsistencies in 
assigning features in the 1000-1200 cm | regions of the IR spectra of humic acids 
either to carbohydrate or to siliceous components. Yonebavashi and coworkers 
(1989) studied forty humic acids from various types of soils by 'H NMR and IR 
spectroscopy. Baes and Bloom (1989) measured spectra of HA and FA from peat 
were obtained using Diffuse Reflectance IR Fourier transform (DRIFT) spectros- 
copy and Fourier Transform-IR (FTIR) spectra using KBr pellets. 

In measuring the IR spectra of HS, it is important to use an extraction procedure 
that does not damage the molecular structure. In performing these extractions, we 
wanted to avoid condensation reactions that may occur. Carbonyl compounds, 
carboxylic compounds, phenols, imines, and aromatic compounds all combine with 
each other and among themselves to form larger molecules. For example, carbonyl 
condensations, or aldol condensations, involve the presence of a carbonyl, an 
enolate, and a strong base. The enolate, a carbonyl compound that acts as an 
electrophile, attacks the carbonyl group under nucleophilic addition to give a p- 
hydroxy ketone or aldehyde also known as an aldol. However, the aldol product may 
further react to give an «,fh-unsaturated carbonyl compound. This is due to the fact 
that the aldol condensation is an equilibrium process. The condensation reactions 
among carbonyl compounds, carboxylic compounds, phenols, imines, and aromatic 
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TABLE |. Correlation table of infrared frequencies for functional groups that could potentially be 
found in humic substances. 


Functional group Wavenumber (cm!) Notes 
C-C] (aliphatic) 600-800 Sharp 
C-Cl (Aryl) 1050-1100 
CH,-Cl (bend) 1200-1300 
P-O 1100-925 
=s 1200-1050 Strong 
C-C 1200 
C-N 1200 Usually strong 
Phenolic OH 1200 Strong 
bending 1410-1310 Strong 
And C-O stretch 
Eo 1300-1250 
C-N 1360-1250 
C-O Alcohols (1220-1040) 
Ethers (1300—1000) 
Aromatic C=C 1650-1450 
N-H 1620-1510 Strong 
-CN 1600 Usually strong 
L=2 1660 
=. 1710 (aldehydes) Very strong 


1715 (ketones) 
1730-1700 (carbox. acids) 
1735 (Esters) 

1655 (amides) 


COOH 1725-1680 Strong 

C=N Just above 2200 Usually strong 
S-H 2600-2550 Weak 

Sp’ C-H 3000-2800 

Acid O-H 3000 Broad 

Sp’ C-H 3100-3000 

Sp C-H 3300 

Alcohol O-H 3300 Broad 

Amine N-H 3300 Broad with spikes 
N-H stretch 

1° 3500 and 3400 Medium 

fam 3500-3310 Medium 
=N-H 3400-3300 Medium 


compounds can take place with the addition of a catalyst to produce a large molecule 
and a smaller molecule, usually water or alcohol. The most common catalysts used 
in these reactions is a strong base such as NaOH, or KOH, however, a strong acid 
or combination of strong acids can also be used for some condensation reactions. 
Because various isolation processes involving HS involve either strong acid or basic 
solutions, we sought a long-term (several weeks) study to identify if any reactions 
took place that would alter the HS structure. 

Past work by Hoekstra and coworkers (1999) has shown that humics can 
become chlorinated under natural conditions. They reported the natural formation of 
4-chlorophenol, 2,4-2,5- and 2,6-dichlorophenol, and 2,4,5-trichlorpenol in four 
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selected areas of a rural Douglas fir forest. Here the humic layer was spiked in situ 
with a solution of Na°’Cl and covered by an enclosure, and incubated for one year. 
Chlorinated phenols were formed naturally from organic matter and inorganic 
chloride presumably by either de novo synthesis or chloroperoxidase-catalyzed 
chlorination. While the conditions are different than those found in the lower 
Suwannee River basin, the study does show that chlorinated aromatics can be 
generated naturally when humics are in the presence of chlorides. 

Recently additional studies on the interaction of humic and fulvic acids with 
chlorine and chloride in soils were published (Casey, 2002; Singer, 1999). These 
studies confirm that organochlorine compounds are present in soil and identify a net 
transfer of chlorine from inorganic to organic forms with weathering. An abundance 
catalytic peroxidase facilitates the chlorination of aromatic compounds by hydrogen 
peroxide in soils. They claim that chlorination in humic materials is widespread and 
may account for organochlorine concentrations found in unpolluted environments. 
Organochlorine compounds attract the attention of environmental scientists because 
of their toxicity and carcinogenic potential in humans. 


MATERIALS AND METHODS—A Mattson Fourier Transform Infrared Spectrometer (FT-IR) was used 
for all measurements. Before using FT-IR to study HS’s from the Suwannee River, we sought to properly 
select conditions for extraction that would minimize any structural damage to the naturally occurring 
organic matter. Infrared spectroscopy has long been used as a tool to qualitatively identify various 
structure aspects (carbonyls, carboxylates, etc.) of HS’s (Yates, 1999; Klavins, 1999; Francioso et al., 
1998). Solutions of deionized ultra filtered (DIUF), 1% HNO3, 20% HNOs:, 0.05 M Na4P,0,, 3 M 
Na>P,07, 0.05 M KOH, and 3 M KOH (100 mL each) were prepared. They were treated with low levels 
of sodium azide to prevent the growth of microorganisms. These solutions were saturated with commer- 
cially available HA (Sigma Aldrich). We then made IR transmission measurements of the seven samples 
at 0, 7, 14, 21, and 28 days. In each case the water was dried and the HS left on the IR transmission card 
was spectroscopically analyzed. We did this to examine what impact, if any, extreme acid or base con- 
ditions would have on HS over a various time periods. 

While our extractions with acids and bases of SRHS were conducted in short time scales (minutes), 
by analyzing longer time scales we hoped to amplify any destructive reactions taking place so the 
insensitive FT-IR could detect these changes. We found substantial structural changes induced by high 
concentration acids or bases and identified appropriate acid (<1% HNO3) and base (<0.05 M KOH) 
concentrations for our work. 

The samples were collected from the sites along the Suwannee River and its tributaries (Table 2). A 
Garmin Handheld Global Positioning System unit was used to mark the latitude and longitude of each 
sample. In each case, approximately 75 cubic centimeters (cm*) of sediment was collected from a site that 
was on the river or spring banks in shallow water. We typically looked for pools of water where there 
appeared to be an accumulation of organic material. The samples were transported back to Valdosta at 
ambient temperature. The material was allowed to settle in a lab (72 °F) for five days. This resulted in 
a sample divided into three layers. The top layer was water and was pipetted off and discarded. The 
second layer was rich in humic material and was removed and placed in a separate container. The bottom 
layer was primarily mineral content (CaCO3, SiO», etc.) and was also discarded. The carbon rich second 


— 


Fic. la,b,c. FT-IR spectrum of humic substances from three different locations on the Suwannee 
River or its tributary. (a) is the IR spectra of HS extracted from sediment taken from the Alapaha River in 
South Georgia. (b) is IR spectra of HS collected on the Suwannee near the Sante Fe River. All freshwater 
samples have similar spectra characteristics. (c) is FT-IR data of HS extract taken in saltwater near the 
Gulf of Mexico, where the appearance of the C-Cl vibrations became apparent. 
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TABLE 2. Description of twenty-six solid (muck!) samples collected throughout the Suwannee river 
basin containing various levels of HS. (s = strong, m = medium, w = weak; b-broad FWHM > 100 
em |, n= narrow FWHM < 100 cm’, c = convoluted with another peak). 


Wavenumber, spectral 


Latitude features, functional 
No. Site Name Description Longitude group assignment 
l Alapaha River Soil sample taken in 30° 42.255 N 599 (m,b) =CH aromatic 
Echols County, six inches of water 83° 01.943 W —- 956 (w,n) C-H rocking 
Georgia from middle of river. 1114 (s,n) C-O stretch 


1681 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 
2924 (m,n) -C-H aliphatic 
3097 (m,c) -CH aromatic 

3324 (s,b) -OH group 


2 Suwannee River, Soil sample taken on 30° 40.987 N 594 (m,b) =CH aromatic 
Fargo, Georgia shore of river in 82° 33.558 W —- 943 (m,n) C-H rock 
shaded area. River 1090 (s,n) C-C stretch 
moving 1-2 knots. 1137 (m,c) C-O stretch 


1669 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

2300 (w,b) -CN 

2920 (m,n) -CH aliphatic 

3090 (m,b) -CH aromatic 

3358 (s,b) -OH group 


3 Suwannee River Taken approximately 30° 40.934 N 493 (w,n) 
Fargo, Georgia 90 meters downriver 82° 33.576 W 582 (m,b) =CH aromatic 
from previous sample. 954 (w,n) C-H rocking 
Sandy bottom with 1126 (s,n) C-O stretch 
organic sediment 1675 (m,b) C=O (ketone, 
sitting on top. etc.), C=C (aromatic) 


2341 (w,b) -CN 

2917 (s,n) -C-H aliphatic 
3068 (m,c) -CH aromatic 
3347 (s,b) -OH group 


4 Withlachochee Low flow but relatively 30° 38.134 N 532 (m,n) 
River, Brooks deep dark water river. 83° 18.703 W 620 (m,b) =CH aromatic 
County, Sample collected 854 (w,n) 
Georgia near bridge upstream 995 (m,n) 
from a paper mill. 1085 (s,n) C-O stretch 
1282 (w,n) 


1658 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

2260 (m,b) -CN 

2491 (m,b) 

2923 (m,n) -C-H aliphatic 

3272 (s,b) -OH group 


5 Withlachochee Down river of paper mill. 30° 35.803 N 592 (m,b) =CH aromatic 
River, near Water darker in color 83° 15.607 W 1103 (m,n) C-O stretch 
Pinetta, Florida than above mill but 1669 (w,n) C=O (ketone, 
minimal green algae ; etc.), C=C (aromatic) 
that was abundant at 2352 (w,b) -CN 
other site. 3361 (s,b) -OH group 


(and CH, NH) 
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TABLE 2. Continued. 


No. Site Name 


6 Suwannee River 
State Park, 
Florida 


7 Suwannee and 
Withlachochee 
River 


8 Withlachochee 
River and 
Suwannee river 


9 Spring at 
Interaction of 
Withlachochee 
and Suwannee 
Rivers 


10 Suwannee River 
south of With. 
River 
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Description 


North (upstream) of 
Withlachochee River 
intersection. Suwannee 
River has carved out a 
limestone channel and 
is moving at low 
(1 knot) rate. Taken 
in 0.2 meters of 
water. 


Intersection of two 
rivers. Took solid 
sample in seven 
feet of water under 
limestone overhang. 


Sample taken on a sand 
bar 0.4 miles north 
of intersection with 
Suwannee River. 


A small spring, 
surrounded 
by a brick wall built 
to hold water, with 
water temperature 
of 72 °F. Soil (muck, 
small twigs, some 
sand, leaves) sample 
taken in 10 feet of 
water in springs 
cave. Sample 
collection site had 
minimal sunlight 
exposure. 


Soil sample taken from 
bank of river 0.65 
miles south of rivers 
intersection. Lots of 
limestone and 
heavily vegetation 
on banks. 


Latitude 
Longitude 


30° 23.581 N 
83° 09.884 W 


30° 23.143 N 
83° 10.266 W 


30° 23.459 N 
83° 10.545 W 


30° 23.206 N 
83° 10.301 W 


30° 22.965 N 
83° 10.791 W 


259 


Wavenumber, spectral 
features, functional 
group assignment 


538 (m,n) 

1081 (m,n) C-O stretch 
1614 (w,b) C=O (ketone, 
etc.), C=C (aromatic) 

2340 (w,b) -CN 

2926 (s,n) -C-H aliphatic 

3318 (s,b) -OH group 
(and CH, NH) 


581 (m,b) =CH aromatic 

940 (m,n) C-H rocking 

1103 (s,n) C-O stretch 

1658 (m,n) C=O (ketone, 
etc.), C=C (aromatic) 

2274 (w,b) -CN 

2920 (s,n) -C-H aliphatic 

3274 (s,b) -OH group 
(and CH, NH) 

592 (m,b) =CH aromatic 

952 (n,w) C-H rocking 

1114 (s,n) C-O stretch 

1658 (w.n) C=O (ketone, 
etc.), C=C (aromatic) 

2341 (w,b) -CN 

3080 (m,n) 

3374 (s,b) -OH group 
(and CH, NH) 

588 (m,b) =CH aromatic 

947 (n,w) C-H rocking 

997 (w,n) 

1114 (s,n) C-O stretch 
1651 (w,n) C=O (ketone, 
etc.), C=C (aromatic) 

2336 (w,b) -CN 
3323 (s,b) -OH group 
(and CH, NH) 


527 (w,n) 

1103 (m,n) C-O stretch 
1636 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 
2926 (m,n) -C-H aliphatic 
3318 (s,b) -OH group 

(and CH, NH) 
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TABLE 2. Continued. 


No. Site Name 


a Peacock Spring 


12 Little Creek 
Spring 


13 Suwannee River 
at intersection 
of Little 
River Spring 


14 Santa Fe River 


5) Suwannee River 
upriver from 


Santa Fe River 
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Description 


Soil samples taken 


from bank at water 
level of spring made 
famous by cave and 
cavern scuba divers. 
Essentially no flow 
from spring and 
entrance quite murky. 


Classic clear water 


springs that empties 
into Suwannee River. 


Solid sample, with lots 


of mineral mixed in, 
was taken in [5 feet 
of water from within 
cave (minimal 
sunlight). Flowing 
with good 

(100 feet) visibility. 
(72 °F) 


Sample taken in about 4 


feet of water about 8 
feet into Suwannee at 
outlet of the Little 
river Spring run. 

The sandy bottom 
had a temperature 
gradient. 


Solid sample taken in 


15 feet of water 

near underwater 
limestone cliff. 
Relatively clear water 
(30 feet visibility) 
with 2 knot current. 


Solid sample taken on 


sand/muck bar near 
the riverbank in 6 
inches of water. 


Latitude 
Longitude 


30° 07.400 N 
83° 07.981 W 


29° 59.802 N 
82° 57.975 W 


29° S97 716IN 
82° 57.995 W 


29° 54.595 N 
82° 51.857 W 


2975328122 N 
822523530 
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Wavenumber, spectral 
features, functional 
group assignment 


538(m,n) 
874 (w,n) 
1081 (s,n) 

1406 (m,n) C=C, CH) 
1614 (m,n) C=O (ketone, 
etc.), C=C (aromatic) 

2374 (w,b) -CN 
3318 (s,b) -OH group 
(and CH, NH) 


581 (m,b) =CH aromatic 

842 (m,n) 

951 (m,n) C-H rocking 

1103 (s,b) C-O stretch 

1372 (w,b) 

1658 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

3318 (s,b) -OH group 
(and CH, NH) 


Sample did not contain 
enough HS to obtain 
measurable IR spectra. 


570 (m,b) =CH aromatic 
842 (w,n) 

951 (w,n) C-H rocking 
1092 (s,n) C-O stretch 
1364 (w,n) 

1669 (w,n) C=O (ketone, 
etc.), C=C (aromatic) 
2934 (m,n) -C-H aliphatic 
3318 (s,b) -OH group 

(and CH, NH) 


479 (m,n) 

581 (m,b) =CH aromatic 

1114 (s,n) C-O stretch 

1680 (m,n) 

2315 (m,b) -CN 

2917 (m,n) -C-H aliphatic 

3318 (s,b) -OH group 
(and CH, NH) 
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TABLE 2. Continued. 


No. Site Name 
16 Suwannee 
River 
t7 Blue Hole 
Spring 
18 Ichnetucknee 


Spring Head 


19 Suwannee 
River 


Second spring located on 


Last set of freshwater 
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Latitude 
Description Longitude 
Solid (muck/mineral) 29° 52.443 N 
sample taken on 82° 52.565 W 
riverbank near cypress 
tree near a limestone 
ledge. 
Spring that gives rise to 29° 58.850 N 
the Ichnetuckee River, 82° 45.505 W 


which empties into 

Santa Fe river. 

Sample collected in 

6 feet of water 10 

feet from springhead. 

Water is 72 °F and 

has >100 feet 

visibility. 

29° 59.077 N 
Ichnetucknee River. 82° 45.693 W 
Mucky sample collected 

20 feet downstream of 

springhead in 1.5 feet 

of water near the 

shoreline. 


29° 35.437 N 
samples taken on bank 82° 56.174 W 
of river approximately 

20 miles upriver of 

Gulf of Mexico near 


Fanning Spring inlet. 
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Wavenumber, spectral 
features, functional 
group assignment 


526 (w,n) 

601(m,b) =CH aromatic 

956 (w,n) C-H rocking 

1103 (s,n) C-O stretch 

1148 (s,n) 

1408 (w,b) C=C, CH> 

1680 (w,b) C=O (ketone, 
etc.), C=C (aromatic) 

2306 (w,n) -CN 

2926 (s,n) -C-H aliphatic 

3318 (s,b) -OH group 
(and CH, NH) 


601 (m,b) 
873 (w,n) 
989 (m,n) 

1431 (w,n) C=C, CH, 
1677 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

2373 (w,b) -CN 
3275 (s,b) -OH group 
(and CH, NH) 


581 (m,b) =CH aromatic 
831 (w,n) 

951 (m,n) C-H rocking 
1114 (s,n) C-O stretch 
1390 (n,w) 

1473 (n,w) C=C, CH, 
1669 (m,n) C=O (ketone, 
etc.), C=C (aromatic) 

2364 (w,b) -CN 

2918 (n,m) -C-H 
aliphatic 

3318 (s,b) -OH group 
(and CH, NH) 


559 (w,n) 
1092 (m,n) C-O stretch 
1451 m,b) C=C, CH, 
1658 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 
2339 (w,b) -CN 
2926 (s,n) -C-H 
aliphatic 
3144 (s,b) -OH group 
(and CH, NH) 


TABLE 2. Continued. 


No. Site Name 


20 Fanning 
Springs 


21 Lower 
Suwannee 
River (near 
town of 
Suwannee, 
Florida)— 
canal 


N 
i) 


Lower 
Suwannee 
River— 
Island 


Dane 
delta-Gulf of 
Mexico 


Solid sample taken near 


Suwannee River —_ Two solid samples taken 


FLORIDA SCIENTIST 


Latitude 
Description Longitude 
Solid sample taken from 29° 35.269 N 
bottom (15 feet depth) 82° 56.140 W 
in clear water (> 100 
feet visibility) spring 
(72 °F). 300-foot 
long spring run ends 
in Suwannee River. 
Solid sample taken in 29° 19.561 N 
1.5 feet of water on 83° 8.343 W 


shoreline of river near 
small town. Murky 
water is brackish 
sometimes. 
Approximately 3 
miles upriver of 

Gulf of Mexico. 


29° 18.525 N 
sea grass bed on 83° 9.125 W 
small island in middle 

of river. Area strongly 

influenced by tides 


(brackish) 


29° 18.283 N 
near channel marker 83° 10.391 W 
11. One taken from 

muck in saw grass, 

and second taken on 

mucky sandy shore on 

shore line of small 

island that has 20 

meter tall oyster pile 

from American 

Indians (age 

500-1500 years 

old). 
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Wavenumber, spectral 
features, functional 
group assignment 


581 (m,b) =CH aromatic 

842 (m,n) 

951 (m,n) C-H rocking 

1105 (s,n) C-O stretch 

1151 (s,n) C-O stretch 

1386 (m,b) C=C, CH, 
CH, 

1669 (m,n) C=O (ketone, 
etc.), C=C (aromatic) 

2380 (w,n) -CN 

2926 (s,n) -C-H aliphatic 

3296 (s,b) -OH group 
(and CH, NH) 


411 (w,n) 

556 (m,n) 

1049 (s,b) C-O stretch 

1415 (w,n) C=C, CH, 

1635 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

2372 (w,b) -CN 

2948 (m,n) -C-H aliphatic 

3296 (s,b) -OH group 
(and CH, NH) 


483 (m,n) 

538 (m,n) 

873 (w,n) 

1049 (s,b) C-O stretch 

1427 (w,n) C=C, CH, 

1625 (m,b) 

2368 (w,b) -CN 

2948 (m,n) -C-H aliphatic 

3296 (s,b) -OH group 
(and CH, NH) 


570 (m,b) =CH aromatic 

633 (m,b) C-Cl (aliphatic) 

879 (w,n) 

985 (m,n) 

1016 (m,n) C-Cl 

1202 (w,n) C-O stretch 

1403 (w,n) C=C, CH> 

1649 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

2278 (w,b) -CN 

2926 (s,n) -C-H aliphatic 

3361 (s,b) -OH group 
(and CH, NH) 
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TABLE 2. Continued. 


Wavenumber, spectral 


Latitude features, functional 
No. Site Name Description Longitude group assignment 
25 Suwannee River Sample taken at low tide 29° 18.383 N 476 (w,n) 
delta-Gulf of on “‘flats” or grass 83° 11.316 W 592 (m,b) =CH aromatic 
Mexico beds in 2 inches of 652 (m,n) C-Cl (aliphatic) 
saltwater near channel 839 (w,n) 
marker 6. 940 (m,n) 
995 (w,n) 


1024 (s,n) C-Cl 

1103 (s,b) C-O stretch 

1452 (w,b) C=C, CH, 

1669 (m,b) C=O (ketone, 
etc.), C=C (aromatic) 

2380 (w,b) -CN 

2980 (m,n) -C-H aliphatic 

3252 (s,b) -OH group 
(and CH, NH) 


layer was passed through a 20-mesh stainless steel screen to remove large particulate matter. We 
collected samples from twenty-five locations along the Suwannee River, its springs and tributaries (see 
Table 2). 

The extraction procedure from purified sample used was as follows. Approximately one gram of the 
cleansed humic material was placed in a container. It was centrifuged and allowed to settle. The aqueous 
layer (top) and mineral content (bottom) were removed, if there was any present. The remaining material 
(middle layer) was washed with 10 mL of 0.05% HNOs3. This was used to dissolve any remaining mineral 
content (i.e., CaCO3) and to remove cations (Mg, Na’, etc.) and anions (SOx, PO; , Cl, etc.) that 
might be present in the HS. The remaining material was washed with 10 mL of deionized water to remove 
any additional ions that might be present. A fraction of the HA and FA present was also removed from 
each sample during the washings. Considering the solubility limit of these species (HA, FA) is in the 10— 
20 ppm range in the aqueous phase, the quantity lost from the sample was quite low (tens of micrograms) 
when compared to the total content in the naturally occurring organic matter. By limiting the quantity of 
water used in the extraction we limited the amount of slightly soluble organic matter separated while 
removing very soluble anions that have vibrations that might interfere with HS spectral features (P-O in 
phosphate/hydrogen phosphate, etc). 

The samples were than analyzed by three FT-IR methods. In the first attempt we used Diffuse 
Reflectance (DRIFT) but found the signal intensity too weak to produce a discernible spectrum. In the 
second approach we extracted the sample with methanol, dried the extracted sample on transmission 
cards, and measured the extract with FT-IR transmission. This approach did not produce a strong enough 
signal to produce a useable infrared transmission spectrum. In the third approach we used DIUF water and 
added enough humic material to make the solution a 200 ppm solution slurry. A few drops of this material 
was placed on an IR transmission card and allowed to dry, in a desiccator, for several days. This last 
approach produced quality spectra that are described in this paper. 


ResuLts—While infrared spectroscopy is not a sensitive technique in terms of 
providing structural information on trace compounds, nor is it a selective, qualitative 
technique in terms of providing us with the molecular composition, it does provide 
us with bulk structural features of our sample. HS are composed of a range of 
naturally occurring organics species derived from the plant and animal decay so 
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when a process such as chlorination has occurred it has impacted a number of 
different molecular species. Some trends observed in analyzing the twenty-six 
samples are outlined above (Table II). 

While broad peaks in the 1600-1700 cm range can be attributed to C=C in 
aromatic structures, there is no spectral evidence in the 1800-2000 cm ' region to 
indicate significant substitution patterns. The convoluted peaks from 2800 to 3500 
cm ‘are intense and broad in all samples. This indicates an abundance and mixture 
of CH aliphatic, CH aromatic, NH», and bonded —OH bonds. 

While carboxylates are widely recognized as the dominant acid group in humic 
and fulvic acid, the relative intensity of this spectral feature (1710 cm‘) has a low 
intensity indicating much of our sample consist of humin. An asymmetric broad 
spectral feature in the 1600-1700 cm” was evident in all samples. This was 
attributed to C=C (aromatic), carboxylates, and carbonyls. The peaks were broad 
and impossible to deconvulte. There was consistent evidence in all samples for C-O 
vibration in the 1030-1220 cm’ regions, assumed to be part of various derivations 
of alcohols, ethers, phenols, etc. There were no strong peaks in the 3500-3600 cm | 
region that would indicate free -OH groups. 

A spectral feature in the 1400-1450 cm ' region was identified in many of the 
samples and was attributed to C-H deformation of CH, or CH groups and aromatic 
C=C. When the peak shifted to 1370-1390 cm", it indicated the presence of a —CH3. 
The spectral features in this region were typically convulted into a single, broad 
(FWHM 75-100 cm‘) profile that was difficult to assign to individual vibrations. 

There was consistent evidence for C=N in the 2300 cm’ region. The peaks 
were typically weak and broad (120 cm | wide) and it was difficult to provide exact 
peak maximum. While they had a relatively low abundance compared to other 
prominent functional groups, the —CN group is attached near a variety of other 
functional groups, meaning we are observing —CN not in a single structure but 
attached to various aliphatic and aromatic environments. We believe this group is 
the product of UV light striking the humics over a period of time. Work by 
atmospheric chemistry has shown that simple organics can react with nitrogen 
containing compounds to produce C=N in the presence of UV light. 

Between 630-680 cm! there are peaks in salt-water samples indicating 
aliphatic chlorides. In the 1240-1220 cm”! region there are spectral features 
indicating aliphatic chlorides in salt-water samples. There were no spectral features 
in the 1800 cm | region from salt-water samples indicating that the chlorination was 
not acid chlorides. There is no evidence for axial or equatorial substituted chlorine in 
a ketone (1750-1730 cm '). A simple analogy can be drawn with the IR spectra 
between 500-1500 cm™' for propane, propanol, chloropropane, and 1-chloro-3- 
propanol. Propane has a single peak in the 1450-1500 cm‘ range. While the exact 
assignments of some of the peaks may be difficult, it can be clearly followed in 
table 3 that a number of spectral peaks in the 1000-600 region appear, strongly, in 
1-chloro, 3-propanol (i.e. 950, 750, 655) and 1-chloropropane (880, 790, 730, 667) 
that are either much weaker or do not appear in the propane and propanol spectra. 
Propane and propanol are simple analogies to the aliphatic characteristics of HS 
found in fresh water. 
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TABLE 3. Infrared Spectral characteristics of propane, propanol, and chlorinated propane (NIST 
Webbook, 2001). 


I 


Species Wavenumber (cm _) 
Propane (C3Hg) 1400-1500 
1-Propanol (C3;H7OH) 1460, 1390, 1215, 1050, 980 
1-chloropropane (C3H7Cl) 1440, 1400, 1300, 1260, 880, 790, 730, 667 
1-chloro, 3-propanol (C3H7CIO) 1505, 1348, 1075, 950, 750, 655, 530 


CoNCLUSIONS—While the humic and fulvic acid fractions of naturally occurring 
organic matter from the Suwannee River at the Fargo, Georgia location has been 
extensively studied by various techniques, we chose to look at humic substances. HS 
comprises a much greater mass percent of naturally occurring organic matter so 
understanding its functional groups and structure can have broad impact in 
understanding its impact on the environment. We also looked at samples from the 
length of the Suwannee River as well as its tributaries and springs. FT-IR is a 
technique that only provides prominent structural features and most of the spectral 
features identified here would be expected, based on past studies, with naturally 
occurring organic matter (C-H, C=C, C-O, etc.). The fact that spectral evidence 
existed for chlorination in two samples indicates that a significant level of 
substitution by the halogen. Since HS are a collection of a large number of relatively 
insoluble organic structures we cannot propose any specific structures or mech- 
anisms at this time. The only structural information that can be derived at this point 
is that substitution took place on an aliphatic structure. We will look into this area in 
the future. 
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EFFECTS OF SUPPLEMENTAL FEEDING ON KEY DEER 
BEHAVIOR IN URBAN AREAS 
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‘Department of Wildlife and Fisheries Sciences, Texas A&M University, TAMU 2258, 
210 Nagle Hall, College Station, TX 77843-2258 
‘National Key Deer Refuge, Big Pine Key Shopping Plaza, Big Pine Key, Florida, 33043 


AssTrAcT: The sociobiology of Florida Key deer (Odocoileus virginianus clavium) appears to have 
changed in recent years. It has been hypothesized that supplemental feeding has increased deer use of 
urban areas. We evaluated and categorized Key deer (n = 132) behavior in urban areas. We found 
feeding was the most frequent (49-54%) behavior observed followed by moving (20-23%), standing (14— 
17%) and bedding (8-12%). In comparing types of feeding behavior, we found feeding on native 
vegetation was more frequently observed (P < 0.001) than other feeding behaviors, including illegal 
feeding by tourists. Illegal feeding during the 6-month tourist season (Oct-Mar) increased (+25%) as 
expected, whereas native feeding behavior decreased (—21%). We propose that illegal feeding is more an 
opportunistic behavior in Key deer rather than the primary reason in use of urban areas. We recommend 
that public outreach and law enforcement efforts discourage illegal feeding of this endangered deer herd 
and that efforts focus on tourists rather than local residents. 


Key Words: Deer behavior, feeding, Florida Key deer, Odocoileus virginanus 
clavium, supplemental feeding, urban development, urban wildlife 


FLORIDA KEY DEER, an endangered sub-species of white-tailed deer in the United 
States, is endemic to the Florida Keys on the southern end of peninsular Florida 
(Hardin et al., 1984). Key deer occupy 20—25 islands within the boundaries of the 
National Key Deer Refuge (NKDR), with approximately 75% of the overall deer 
population on Big Pine and No Name keys (Lopez, 2001, Lopez et al., 2003). Urban 
development in the Keys has been a concern in the recovery and management of 
Key deer for the last 3 decades (Klimstra et al., 1974). Since 1970, for example, the 
human population has increased nearly 10-fold on Big Pine and No Name Keys; 
in 2000, the resident population was estimated at 5,000 people with an influx of 
1,000—1,500 tourists during winter (Folk, 1991; Lopez et al., 2003). As a result, the 
sociobiology of Key deer appears to have changed overtime, whereas deer now 
exhibit use of urban areas, larger group sizes, and a reduction in deer movements 
(Folk and Klimstra, 1991; USFWS, 1999). Supplemental feeding of Key deer is 
prohibited by state (F.A.C. 39-27.002) and federal (16 U.S.C. 1531) laws (USFWS, 
1999). Folk and Klimstra (1991) hypothesized that supplemental feeding had 
increased deer use of urban areas in recent years. Reported impacts due to illegal 
feeding include a “‘taming” or urbanization of Key deer (Folk and Klimstra, 1991) 
and an increase in disease transmission (Nettles et al., 2002). However, the role of 
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illegal feeding on Key deer behavior in urban areas is poorly understood. 
Furthermore, the effect of tourism fluctuations on Key deer behavior in urban areas 
also is largely unknown. The objective of our study was to evaluate Key deer 
behavior in urban areas. Specifically, we attempted to (1) determine the type and 
relative frequency of Key deer behavior in urban areas and between seasons (tourist 
vs. non-tourist), and (2) evaluate the impacts of illegal feeding on Key deer feeding 
behaviors in urban areas. 


Stuby AREA—The Florida Keys are a chain of small islands approximately 200 km from peninsular 
Florida. Big Pine (2,548 ha) and No Name (461 ha) Keys are within the boundaries of the NKDR, Monroe 
County, and supports most of the Key deer population (Lopez, 2001). All data presented in this paper were 
collected on these 2 islands. Typically, island areas near sea level (maritime zones) are comprised of red 
mangrove (Rhisophora mangle), black mangrove (Avicennia germinans), white mangrove (Laguncularia 
racemosa), and buttonwood (Conocarpus erecta) forests. With increasing elevation, maritime zones 
transition into hardwood [e.g., Gumbo limbo (Bursera simaruba), Jamaican dogwood (Piscidia 
piscipulaa)] and pineland [e.g., slash pine (Pinus elliottii, saw palmetto (Serenoa repens)| upland forests 
intolerant of salt water (Dickson, 1955). Approximately 24% of the native vegetation has been converted for 
residential or commercial uses since 1955 (Lopez, 2001). 


MetHops—Key deer (n = 206, males = 90, females = 116) were live-trapped and marked between 
February 1998-January 1999 on Big Pine and No Name Keys. Deer were captured using portable drive 
nets (Silvy et al., 1975), drop nets (Lopez et al., 1998), and hand capture (Silvy, 1975). Sex, age, capture 
location, body weight, radio frequency (if applicable), and body condition were recorded for each deer 
prior to release. Captured deer were marked in a variety of ways depending on sex and age. Methods 
included plastic neck collars (8-cm wide, primarily females of all age-classes), leather antler collars (0.25- 
cm wide, yearling and adult males only), or elastic expandable neck collars (3-cm wide, primarily male 
fawns/yearlings, Lopez et al., 2003). Approximately 72% of deer (149/206) were radio-collared with 
a battery-powered mortality-sensitive radio transmitter (150—152 MHz, 100-110 g for plastic neck collars, 
10-20 g for antler transmitters and elastic collars, Advanced Telemetry Systems, Isanti, MN). 

We observed Key deer behavior weekly from February 1998—January 1999 during daylight hours. 
Deer observations were not recorded during evening hours because of concerns with the privacy of 
subdivision residents. Four subdivisions were selected on Big Pine (Port Pine, 46 ha; Koehn, 29 ha; Eden 
Pines, 23 ha) and No Name (Bahia Shores, 14 ha) keys to record deer observations. Subdivisions were 
selected based on high deer densities and the occupancy of marked deer. Human density, pattern of 
development, and deer density were similar among these 4 subdivisions. In recording deer observations, 
we first randomly selected a subdivision. Once a subdivision was selected, a randomly-generated list of 
marked deer was used to identify a focal individual. If the first deer on the list was not found after an 
extensive search (i.e., driving subdivision, use of radiotelemetry), searching of subsequent deer on the list 
was conducted until a focal individual was identified. Once found, marked deer were observed at 
a distance (20-75 m) with the aid of 8 X 35 binoculars from a stationary vehicle. Observations were 
recorded while the marked deer was in sight and in the subdivision. 

We recorded deer behavior every 15 sec to generate a relative frequency (i.e., “time budget’, 
Altmann, 1974; Grier and Burk, 1992) estimate of Key deer behavior in urban areas. We classified deer 
behavior into 4 general categories: standing, moving, feeding, and bedding (Hardin, 1974). Feeding 
behavior was classified into 4 sub-categories: native, planted, illegal feeding, and garbage. We defined 
native feeding as deer eating native vegetation in vacant lots or roadsides. We defined planted feeding as 
deer eating landscape or ornamental plants in urban areas. Illegal feeding behavior was defined as deer 
eating from food given by or put out by humans. Finally, garbage feeding was defined as deer eating from 
garbage cans including vegetation clippings. 

We divided the calendar year into 2 seasons: (1) tourist season (6 months, Oct—Mar) and (2) non- 
tourist or local season (6 months, Apr-Sep). We tested (P < 0.05) for differences in deer behaviors among 
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Fic. 1. Key deer behavior (%) in urban areas between males and females, 1998-1999. 


categories and seasons using Chi-square tests (Ott, 1993). We pooled data between age-classes due to 
small sample sizes, and when appropriate (White and Garrott, 1990). 


RESULTS—We collected observations on 132 deer in our study (female = 87, 
male = 45; 14,400 interval point values, 60 hrs, 10 min). Observation times ranged 
from 2 min to 2 hrs, 13 min. Mean observation time was approximately 20 min. 
Only 4 observational periods were less than 5 minutes and were not removed from 
our analysis because (1) of the potential loss of important behavioral information 
(e.g., deer moving out of a neighborhood), and (2) these 4 deer observations were 
a relatively small proportion of total sample size. 

We found differences (X” = 173.90, d.f. =3, P < 0.001, Figs. 1—2) in standing, 
moving, bedding, and feeding behaviors between males and females. Feeding was 
the most frequent deer behavior (49-54%, ranges represent differences between 
males and females) observed followed by moving (20-23%), standing (14-17%) 
and bedding (8-12%, Fig. 1). We also found differences (X* = 32.96, d.f. =3, P < 
0.001, Fig. 3) in all feeding behaviors between males and females. We found feeding 
on native vegetation (43-64%, ranges represent differences between males and 
females) was more frequently observed than other feeding behaviors (planted 15— 
19%, illegal feeding 4-29%, and garbage 12-13%, Fig. 3). Individual Chi-square 
tests were used in comparing feeding behaviors by season (sexes were combined for 
this comparison); we found differences Oe = 4521 22,00" dt — le P = 0.001) tor 
native, planted, and illegal feeding, however, we found no difference OG = 0:64, 
d.f. = 1, P=0.423, Fig. 4) in the frequency of observed garbage feeding behavior. 
Illegal feeding during the tourist season increased (+25%) while native (—21%) and 
planted (4%) feeding behaviors decreased. 
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Fic. 2. Key deer behavior (%) in urban areas between the tourist season (Oct-Mar) and non-tourist 
season (Apr-Sep), 1998-1999. 


Discussion—We found approximately 50% of observed Key deer behavior 
in urban areas was associated with feeding. Other studies have documented the use 
and importance of urban areas to Key deer (Hardin, 1974). Some of the reported 
urban-related benefits to deer include an increase of food resources due to dis- 
turbances and planting of ornamental plants (Hardin, 1974; Klimstra et al., 1974; 
Hardin et al., 1984, Lopez, 2001). For example, Hardin (1974) reported feeding (his 
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Fic. 3. Feeding behavior (%) in urban areas between male and female Key deer, 1998-1999. 
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Fic. 4. Feeding behavior (%) for Key deer in urban areas between the tourist season (Oct-Mar) and 
non-tourist season (Apr—Sep), 1998-1999. 


study, 38%; our study, 49-54%) was a frequent behavior; however, it should be 
noted that results presented by Hardin (1974) included Key deer observations in 
natural areas such as hammocks and pinelands, as well as developed areas. Thus, 
assuming that feeding on urban vegetation is greater than other vegetation types 
(e.g., buttonwoods), the lower frequency of feeding reported by Hardin (1974) could 
be attributed to the data pooling in his study. 

In comparing feeding behavior between males and females, we found feeding 
on native vegetation (43-64%) was more frequently observed than other feeding 
behaviors (Fig. 3). Study results contradict what is often cited as a primary reason 
for Key deer use of urban areas—illegal or supplemental feeding (Folk, 1991; Folk 
and Klimstra, 1991). Our data suggest illegal feeding as reported by Folk (1991) and 
Folk and Klimstra (1991) is not the primary reason for use of urban areas by Key 
deer. However, we did find illegal feeding during the tourist season increased while 
native feeding behavior decreased. The opposite was found during the non-tourist 
season (Fig. 4). We propose illegal feeding is more an opportunistic behavior in Key 
deer rather than the primary reason in use of urban areas. Observed behavioral 
differences between feeding on native vegetation versus illegal feeding between 
seasons suggests a need to educate tourists rather than local residents about the 
negative impacts of supplemental feeding on this endangered deer population. 
Supplemental feeding should continue to be discouraged; however, public outreach 
and law enforcement efforts may want to concentrate their efforts with tourists rather 
than local residents. 
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ABSTRACT: A stage-structured matrix model for the endangered American Crocodile (Crocodylus 
acutus) in Florida was constructed using published data and biologically reasonable estimates. The 
baseline version of the model yields a population growth rate Aq of 1.006, consistent with observations 
that the C. acutus population is slowly increasing in Florida. Elasticity or the proportional sensitivity of 
iq with respect to changes in annual survivorship indicates that annual survivorship of sub-adults (1 m— 
2.25 m) has the greatest proportional effect on population growth. Annual survivorships can be affected 
by management and are, therefore, the parameters of most interest in this elasticity analysis. Simulated 
management scenarios further demonstrate the degree to which sub-adult annual survivorship can affect 
population growth. The finding that sub-adult annual survivorship has the greatest proportional effect on 
dq suggests that management efforts should be concentrated toward this stage of the American crocodile. 


Key Words: Matrix model, American crocodile, Crocodylus acutus, conservation, 
elasticity 


THE endangered American crocodile (Crocodylus acutus) is a broadly ranging 
species in much of tropical America, including a small and isolated population 
in southern Florida. This population has probably never been larger than 2,000 
individuals since the beginning of the 20th century (Ogden, 1978). The population 
has declined due to habitat loss and human disturbance but is now slowly recovering 
(Ogden, 1978; USFWS, 1984, 1999). Current estimates range from 100-400 non- 
hatchlings (Ogden, 1978) to 500 individuals overall, of which 220 + 78 are non- 
hatchlings (Kushlan and Mazzotti, 1989b). 

There are a variety of hypotheses regarding the key factors affecting the recovery 
of C. acutus in Florida. Among those published are human disturbance (Kushlan and 
Mazzotti, 1989a; Moler, 1991), hurricanes (Ogden, 1978), interactions with alligators 
(Alligator mississippiensis; Moore, 1953), nest site availability, and desiccation and 
flooding of nests (Mazzotti et al., 1988; Mazzotti, 1989). A number of studies have 
also shown that salinity levels affect the survivorship of hatchling and juvenile 
C. acutus (Evans and Ellis, 1977; Ellis, 1981; Dunson, 1982; Mazzotti and Dunson, 
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1984; Dunson and Mazzotti, 1989). However, the relative importance of hatchling 
and juvenile survivorship on population growth is unknown. Of all these effects, the 
most important may be human disturbance, which can manifest itself in many ways. 
While a number of individuals are killed on highways (Kushlan, 1988; Moler, 1991), 
a stable and growing population of C. acutus occurs in the disturbed cooling canals of 
the Florida Power and Light Company’s Turkey Point nuclear power plant (Gaby 
et al., 1985). The relative importance of life history stages to population growth of 
C. acutus 1s poorly understood and needs to be investigated. 

Modeling can relate the effects on individual life stages to the whole population, 
and thus assist in the identification of critical life history stages that can be targeted 
by management. These life stage specific effects can be investigated using stage 
structured models, such as matrix projection models. Matrix population projection 
models based on age were first introduced to biology by Bernardelli (1941), in- 
dependently generalized by Lewis (1942) and Leslie (1945), and later generalized to 
stage-based models by Lefkovitch (1965). Stage-based population projection matrix 
models (Lefkovitch matrices) are particularly well-suited to long-lived organisms, 
e.g. loggerhead sea turtles (Crouse et al., 1987; Crowder et al., 1994), and killer 
whales (Brault and Caswell, 1993), including crocodiles, because age-specific demo- 
graphic data are difficult to obtain and because survival and reproduction depend 
more upon size/stage than upon age. 

Population modeling of crocodilians has concentrated primarily on simulating 
harvest strategies, not on the relative impacts of management strategies. Examples 
include harvest models of A. mississippiensis (Nichols et al., 1976; Rice, 1996), 
C. niloticus (Blomberg et al., 1982), and C. johnstoni (Smith and Webb, 1985). These 
models used large data sets and made fairly accurate predictions of population size. In 
one case (C. johnstoni simulation model, Smith and Webb, 1985) a sensitivity analysis 
was conducted on the input parameters. Their results showed that predicted population 
size was most sensitive to the measured error in egg survivorship, but this parameter 
also had the largest standard error of all the input parameters. They concluded that 
a better population size estimate would be gained by increased attention to the mea- 
surement of egg survivorship, presumably by lowering the measurement error 
associated with it, but they warned that using their results to assign relative importance 
to input parameters is ““misleading”’. A proportional sensitivity analysis (elasticity) 
based on hypothesized true parameter values and not on measurement error would be 
needed to draw such a conclusion. 

In this paper, I use a variety of published data to develop a stage-based pop- 
ulation projection matrix model for C. acutus in Florida. Three primary areas are the 
current focus of C. acutus activity in Florida; the southern portion of Everglades 
National Park, Key Largo, and the Florida Power and Light Company’s Turkey 
Point power plant. Data from published studies in all 3 areas were combined to 
parameterize the model. Using a prospective analysis (Horvitz et al., 1997), I 
addressed the question: which of the model parameters would have the largest effect 
on population growth? Elasticity analysis identifies the parameters having the great- 
est proportional effect on population growth and simulated management scenarios 
demonstrate these effects on C. acutus populations. 
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Fic. |. Life cycle graph of C. acutus. P;’s are the probabilities of remaining in stage classes, G;’s 
are the transition probabilities of moving to the next stage classes and F;’s are the number of hatchlings 
produced per female of stage 1 per year. 


MetTHOoDs—The model—A stage-based matrix model for C. acutus was developed following the 
methods of Caswell (1989). Stage categories based on age, size (total length), and reproductive capacity 
were chosen for biological reasons in a similar way that has been done for plants (Horvitz and Schemske, 
1995) and turtles (Crouse et al., 1987; Crowder et al., 1994; Heppell et al., 1996). Although approximate 
size and stage classes for C. acutus are referred to in the published literature, no general convention exists. 
Mazzotti (1983) and Kushlan and Mazzotti (1989b) referred to hatchlings (<0.5 m), subadults (0.5—2.5 m) 
and adults (>2.5 m). Gaby and co-workers (1985) used size classes of juvenile (60—120 cm), sub-adult 
(121-183 cm) and adult (>183 cm). From these categories and data from Moler (1991) I defined the 
following five stages for C. acutus : 1) “hatchlings” (<0.5 m, <1 year), 2) “juveniles” (0.5 m—1.0 m, | to 
<3 years), 3) “sub-adults” (>1.0 m—2.25 m, 3 to <10 years), 4) “‘young-adults”’ (>2.25 m—2.5 m, 10 to 
<12 years), and 5) “adults” (>2.5 m, 12 to 25 years). The hatchling stage defined here is consistent with 
Mazzotti (1983) and Kushlan and Mazzotti (1989b), but I divided their sub-adult stage into three stages. 
This division minimizes the probability of violating the assumption that all individuals within a stage have 
the same transition probabilities (Moloney, 1986). In Moler’s (1991) 10-year study, C. acutus exhibits 
a change in annual survivorship at about 3 years of age, corresponding to a length of approximately 1.0 m. 
This split at the 1.0 m size results in two stages, juvenile and sub-adult. The “‘young-adult”’ stage was added 
because it has been observed that small (2.25 m—2.5 m), female C. acutus have reduced clutch size (Moore, 
1953; Mazzotti, 1989). Finally I adopted the adult stage of Mazzotti (1983), Kushlan and Mazzotti (1989b) 
and others. 

Figure | shows the five stage model life cycle graph for C. acutus where P; is the probability of 
surviving and remaining in a stage 1, G; is the probability of surviving and growing into the next stage, and 
F; is the fecundity of the stage 1. Writing this as a matrix yields the basic form of the projection model for 
C. acutus (Eqn. 1): 


OgnOr HOWE. iE 
Gremio 0, 10 

A= lO Gi Ox (1) 
Onan G =P, 1-0 
One RoltonG, =P. 


The diagonal (P;) and subdiagonal (G;) elements of A were calculated as a function of annual 
survivorship, p;, and stage duration, d;, (see Parameterization below and Table | for C. acutus estimates). 
P; is calculated (Egn. 2): 
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TABLE |. Parameter sets for version | (the best estimate), version 2 (the minimum) and version 3 
(the maximum). Units are as follows: size is total length in meters, d; is in years, p; has units of per year, 
and F; are in units of hatchlings per year. 


Stage di Pi F 
description Size range (stage duration) (survivorship) (fecundity) 


a) Version | parameters. 


1) hatchlings <.5m l 0.2040 0 
2) juveniles (0.5 m—1.0 m) 2 0.6502 0 
3) sub-adults (>1 m—2.25 m) oi) 0.7862 0 
4) young adults (2.25 m-—2.5 m) 2 0.90 es) 
5) adults = 2.5m 25 0.90 163% 
b) Version 2 parameters. 
1) hatchlings << Sail l 0.0680 0 
2) juveniles (0.5 m—1.0 m) 3 0.3640 0 
3) sub-adults (2 1im=22250m) 9 0.7076 0 
4) young adults (>2.25 m-—2.5 m) 5 0.80 O92 
5) adults 22.) m 5 0.80 1.45 
c) Version 3 parameters. 
1) hatchlings <5 ml it 0.4300 0 
2) juveniles (O' m=—120'm) l 0.7152 0 
3) sub-adults (>1 m—2.25 m) 6 0.8648 0 
4) young adults (>2.25 m-—2.5 m) 0.99 15.76 
5) adults > 2.9m 40 0.99 15.44 
and G; (Eqn. 3): 
1 
pi'(1 — pi) 
=P (3) 
— Pi i 


(Crouse et al., 1987). The population rate of increase Ag (the dominant eigenvalue of A), the stable stage 
distribution (the right eigenvector of A associated with Ag), and the reproductive value vector (the left 
eigenvector of A associated with Ag) were calculated using EISPACK, a public domain subroutine. 
Proportional sensitivities or elasticities measure the relative effects of proportional changes in 
various parameters on Ag. Elasticity analysis is useful since it allows a comparison of the effects of various 
parameters that may not be on the same scale. The elasticities of 44 for each parameterization of the model 
to their respective matrix elements (P;, G;, and F;) were calculated using the computational form found in 
Crouse and co-workers (1987). Since the matrix elements are functions of annual survivorship (p;) and 
stage duration (d;), the elasticities of Ag with respect to these parameters can be calculated (Eqn. 4): 


ha(x + 0.01x) — Ag(x — 0.01x) 


Ee 
0.02A4 


(4) 
(Crowder et al., 1994), where Ay(x + 0.014) is the population growth rate Ay calculated when parameter x 
(either p; or d;) is increased or decreased 1%. This form is an approximation, not an exact analytical form, 
as is the equation for elasticity to matrix elements given in Crouse and co-workers (1987). It was checked 
against a computational form (Wolff, 1997) and yielded nearly identical results that are accurate to 
approximately 4 decimal places. 


Parameterization—General considerations—\ adopted the “‘no holds barred methods of parameter 
estimation” (Caswell, 1989), using a mixture of methods to acquire parameter estimates. Most parameter 
estimates were based on a variety of published values, except annual survivorship (p;) and stage duration 
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(d;) of young adults and adults (see Stage duration and Annual survivorship sections below). Three versions 
of the model were developed using three sets of parameters that span the range of likely values. The 
baseline version | of the model was based on the means of published values for C. acutus in Florida, and 
my best estimates (Table 1a). In addition, minimum and maximum sets of parameters were created either by 
using the extremes of the reported range of values, or by adding or subtracting 10% of the parameter values 
used in version |. This was done for all parameters except the stage duration (d;) of hatchlings, juveniles, 
and sub-adults (see Stage duration and Annual survivorship sections below). The combination of the 
minimum and maximum parameter sets that yield the smallest and largest 4's, were labeled version 2, the 
‘minimum’, and version 3, the ‘maximum’, projection matrices. Table | includes all 3 parameter sets I used 
for analysis. Similar methods of combining means, point estimates and informed guesses have been used 
for other demographic models where few data were available, e.g., a Ridley sea turtle matrix model 
(Heppell et al., 1996), a desert tortoise model (Doak et al., 1994) and a C. johnstoni simulation model 
(Smith and Webb, 1985). 


Fecundity, F;—-The census interval for the model is one year, with censusing in late August 
immediately after hatching. During a census interval, a female will mate, build a nest, and lay eggs, which 
hatch just prior to the next census. Stage 5 (“‘adult’’) fecundity (F; in Figure 1) was estimated using Kushlan 
and Mazzotti’s (1989b) summary of previously published (up to 1982) and unpublished data on the nesting 
biology of C. acutus in Florida. Using their estimates of a mean clutch size of 38 eggs (range 15-56), 
a probability that a nest is successful of 0.783, the proportion of eggs that hatch given a successful nest of 
0.575, and adult female survivorship during the interval, I estimated that 15.25 female and male hatchlings 
were produced per nesting female. Although Deeming and Ferguson (1989) claimed that female-biased sex 
ratios are likely in most crocodilian species, Thorbjarnarson (1997), using published data on 11 species of 
crocodilians, showed that a female-biased sex ratio is not supported. I assumed a |:1 sex ratio for C. acutus 
in Florida (Kushlan and Mazzotti, 1989b; Moler, 1991), and since the model only follows females I divided 
the number of hatchlings per breeding female by two. This gave the estimated female hatchlings produced 
per reproductive female as F; = 7.63 for C. acutus in Florida (Table 1a). Finally, while it is possible that 
female C. acutus do not nest every year, currently there are no data to support an estimate of the proportion 
of females that do nest and, as will be apparent, the range of values for fecundity that I use in the mintmum 
and maximum versions of the model would likely take this into account. 

For the minimum and maximum estimates of F; I attempted to span the likely range of feasible 
fecundities for C. acutus. Using the ranges in Kushlan and Mazzotti (1989b), and including Ogden’s (1978) 
estimate of 0.65 for the probability of a successful nest and female survivorship during the interval, 
a minimum fecundity estimate of 1.45 and a maximum of 15.44 were calculated (Table 1). Published 
data for C. acutus from other parts of its range report: 27 eggs in a single nest in Mexico (Aguilar and 
Casas-Andrew, 1991), a range of 19 to 41 eggs per nest in Colombia (Abadia, 1996), a mean of 22 eggs per 
nest (range = 9—36) in the Dominican Republic (Schubert et al., 1996), and a total fecundity value of 16.04 
male and female hatchlings per nesting female in Haiti (Thorbjarnarson, 1988). These observations fall 
within the ranges that I used for C. acutus in Florida. In addition, a mean of 36.4 eggs per nest (n=21, SD= 
13.93, range = 15—75) was reported by Greer (1975) for C. acutus but the locations of the nests were not 
reported. 

The fecundity estimates (F4, Fig. 1) for stage 4 (young breeders) were calculated similarly. Mazzotti 
(1989) observed a clutch of 22 eggs laid by a 2.28m female and Moore (1953) observed two captive female 
C. acutus of approximately 2.28 m laying clutches of 21 and 22 eggs. Given these limited data, I used 
a clutch size of 22 eggs per stage 4 female and the probability of any egg hatching of 0.450 (see above). 
After accounting for female survivorship during the interval, and growth into the large stage, this yielded 
a value of F; = 5.99 female hatchlings produced per stage 4 female per year for the version | model. For the 
minimum and maximum parameters of F, I reduced the appropriate estimates of the number of hatchlings 
produced per nest of stage 5 females by 50% and then adjusted these by survival and growth of stage 4 
(Table 1). 


Stage duration, d;—The approximate durations of each stage (d;’s of Table 1) were estimated for 
stages 1-3 using Moler’s (1991) Figure 3. I approximated the age of first reproduction (at 2.25 m) to be 10 
years based on Moler’s (1991) Figure 3 and the stage durations of stages 1-3. Age of first reproduction 
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was estimated to be 13 years by both Lebuff (1957) and Ogden (1978) but Moler’s (1991) data suggest 
that C. acutus can reach reproductive size (2.25 m) at age 9 years. Using this possible range of 9 to 13 
years for the age of first reproduction, I set the stage durations for stages 2 and 3 in the minimum and 
maximum versions of the model (Table 1). For stages 2, 3 and 4 the minimum version has the longest 
stage durations and the maximum version has the shortest. This is because, for example, the greater the 
amount of time individuals remain in early life stages the less likely they are to survive to breed at full 
potential and the later they will start reproduction, thus reducing the population growth rate (Aq). 

Longevity in crocodilians is notoriously hard to estimate, so I used a range that is reasonable. For 
version | of the model, the duration of stages 4 and 5 were arbitrarily chosen to be 2 and 25 years, 
respectively. For version 2, the minimum parameter set, I allowed stage 4 to have a duration of 5 years and 
stage 5 to have a duration of 5 years. For version 3, the maximum parameter set, stage 4 has a duration of 
1 year and stage 5 is 40 years. The maximum possible life duration is 37 years in version 1, 23 years in 
version 2, and 49 years in version 3. Although there are no data on the longevity of adult C. acutus in the 
wild, there are accounts of other species of crocodilians living 40-56 years or more in captivity (Levy 
1991), and thus 49 years seems to be reasonable for C. acutus in Florida. 


Annual survivorship, p;—Age-specific annual survivorship data, using minimum number known 
alive, is available for C. acutus from a 10 year mark-recapture study on Key Largo (Moler 1991). I converted 
Moler’s (1991) age-specific annual survivorship data to stage-specific annual survivorship, by assuming that 
annual survivorship within a stage is independent of age. This resulted in estimates for the annual 
survivorship of stages | to 3. Because I am not aware of any published data for annual survivorship of stages 
4 and 5, I chose the midpoint of Moler’s (1991) estimate that C. acutus older than 5 years have annual 
survivorships between 0.80 and 1.00 (0.90). 

To arrive at minimum and maximum estimates of annual survivorships I used other authors opinions 
and extreme estimates, when available. For hatchling annual survivorship I used the extremes of the range 
of 0.068—0.429 (Moler, 1991). The low value is consistent with both Brandt and co-workers (1995), who 
reported that over a 9 year period the mean hatchling annual survivorship was 0.085 at the Turkey Point 
site, and Mclvor and co-workers (1994) who report that annual survivorship of hatchlings in Florida Bay is 
0.10. For stages 2 and 3 I arbitrarily added 10% or subtracted 10% from the version | estimates for annual 
survivorship except the minimum for stage 2. The minimum for stage 2 (0.364, Table 1) was calculated 
using the data of Brandt and co-workers (1995: Table 1). For adults (stages 4 and 5) I set the range = (0.80— 
0.99) (modified from Moler, 1991) for annual survivorship for the mintmum and maximum estimates of 
version 2 and version 3 of the model respectively (see Table 1). 


Projection matrices—The population projection matrices for all three versions of the models were 
calculated from the appropriate parameter sets of Table |. The best estimate, version 1, projection matrix 
is presented in Table 2. The population projection matrices for versions 2 and 3 are readily calculated from 
the parameter sets of Tables 1b and Ic, respectively. 


Management simulations—Four simulated management scenarios (simulations A—C) were created 
as examples of the types of questions that can be addressed by this model and to provide an example of 
what the elasticity results mean. All simulations began with a total population of 500 individuals 
distributed according to the stable age distribution produced by version | of the model and were projected 
for 20 years by repeated matrix multiplication. Simulations A and B represented a hatch and release 
program by increasing fecundity 21% and maintaining that level of fecundity for 2 and 20 years 
respectively. Simulation C involved increasing sub-adult (stage 3) survivorship from 78% to 80.9% for all 
20 years of the simulation, to simulate the effects of reducing road mortality. The fecundity increase in 
simulation B was chosen to produce a population growth rate close to that of simulation C. 


ResuLts—Eigenvalue analysis and elasticities—Each of the three projection 
matrices (version | the ‘best estimate’, version 2 the ‘minimum’ and version 3 the 
‘maximum’) yielded a dominant eigenvalue (Ag) (population growth rate) that 
ranged from 0.726 for version 2 to 1.377 for version 3. A model with this much 
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TABLE 2. The population projection matrix for version 1 of the model. 


Stage 
Stage l 2 3 4 5 
1 0 0 0 Des/y) 7E68 
2 0.204 0.394 0 0 0 
3 0 0.256 0.738 0 0 
4 0 0 0.049 0.474 0 
5 0 0 0 0.426 0.892 


uncertainty is not useful for making quantitatively accurate predictions for the 
population growth rate. However, it can be used to examine the proportional effects 
of the model parameters on population growth, or elasticity of Ay to changes in the 
parameters. 

The general pattern of elasticities for all matrix elements (P;, G;, and F;) in all 
three versions of the model were similar, with one minor difference; thus, I present 
the elasticities of Aq only for the matrix elements of version | (Flg. 2). For all 
versions the elasticities to P3 and Ps are larger than all other elasticities to matrix 
elements except for the elasticity to Py in version 2. This difference is due to the 
effect of the stage duration of stage 4 (dy) on the elasticities of Ag to P4. In version 1, 
where d4 = 2, it is expected that 4, would be more sensitive to changes in P, than 
when d, = | in version 3, since P4 = O and elasticity is zero. Likewise, the elasticity 
of A, is still higher in version 2, where dy =5. This result disappears if I set d4 = 3 in 
version 2. These results showed that relatively small changes in P3 and Ps; cause 
much greater changes in the population growth of C. acutus than proportionally 
similar changes to any of the other matrix entries. 

Since the matrix elements were functions of the lower level parameters, d; and 
p;. elasticities of Ag to these parameters may be more useful for assessment of 
management decisions on population growth (Crouse et al., 1987). Elasticities of Ay to 
annual survivorship (p;) are presented in Flgure 3 for all three versions. The general 
pattern of elasticities of Ag to p;'s were similar for all three versions, and the qualita- 
tive results were the same. Sub-adult annual survivorship (p3) had the largest 
elasticity of Ag in all three versions. Adult annual survivorship also showed high 
elasticities of Ag. 

Elasticities of Ag to stage duration were all small and usually negative. All 
elasticities to stage duration for all 3 versions fell within the range of —0.194 to 
0.074, which are comparable to the values found by Crouse and co-workers (1987) 
in their elasticity analysis on loggerhead sea turtles and were small relative to the 
elasticities of annual survivorship for stages 3 and 5 (Fig. 3). 


Results of simulations—Version 1 of the model was projected for 20 years to 
provide a baseline projection to compare to simulations A—-C. Flgure 4 shows the 
results of the 20-year projections in terms of the non-hatchling population. Version 1| 
has a Ag = 1.006 which indicates slow population growth; the simulated population 
increased from 220 to 248 non-hatchlings in 20 years. Figure 4 also includes the 
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Fic. 2. Elasticities of Ag with respect to the matrix elements for version 1. Versions 2 and 3 have 
a similar qualitative pattern of elasticities, with the largest elasticities due to P;’s for stages 3 and 5, except 
for the addition of a large elasticity due to Py in version 2 (see text for details). 


results of the four simulated management scenarios (A—C). Of these, increasing sub- 
adult annual survivorship to 80.9% (simulation C) yielded the highest population 
growth rate, Ag = 1.0192, while increasing fecundity 21% over all 20 years 
(simulation B) yielded a similar Ag (1.0188). At the end of the 20 year projection 
increasing fecundity 21% (B) resulted in a greater total non-hatchling population 
(see Fig. 4) but it would eventually be surpassed by simulation C if the projection 
length was increased sufficiently. The simulated 2 year hatch and release program 
(A) resulted in a short term increase in the non-hatchling population but soon returns 
to nearly baseline levels resulting in only a 2.8% increase over baseline (version 1) 
in the simulated population by year 20. 


DiscussiION—Elasticity analyses and simulations using the best estimate version 
of the model of the American crocodile (C. acutus) strongly suggest that annual 
survivorship of sub-adults (individuals 1 m—2.25 m) has the greatest effect on 
population growth rate. The same qualitative result was found in all versions of the 
model (Fig. 2). This conclusion is supported by similar models and empirical data 
that show the greatest effect on population growth rate of sub-adult or young adult 
survivorship for other large, long-lived species such as sea turtles (Crouse et al., 
1987; Crowder et al., 1994; Heppell et al., 1996), desert tortoises (Doak et al., 1994), 
killer whales (Brault and Caswell, 1993), and long-lived birds, e.g. the spotted owl 
(Lande, 1988) and Hawaiian dark-rumped petrel (Simons, 1984). Because annual 
survivorship is a parameter that management can affect, it is more appropriate to 
examine the elasticity of Ag to this parameter than to matrix elements which are 
functions of both annual survivorship and stage duration, such as the elasticity of 
hq to the probabilities of sub-adult and adult stasis (P;, Fig. 2). In any case, for 
C. acutus in Florida, management practices that may increase annual survivorship 
of sub-adults (e.g. reducing road mortality) would probably also increase annual 
survivorship for all stages >1 m total length, which would increase Ps as well. 
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Fic. 3. Elasticities of Ag with respect to annual survivorship (p;) for all 3 versions. Elasticities are 
greatest for stage 3 sub-adults for all 3 versions. 


Making quantitative predictions from models is risky (Crowder et al., 1994), and 
this is certainly true of this model. The estimates for adult annual survivorship are only 
rough estimates, and most of the remaining parameter estimates were derived from 
only a handful of studies. I assumed a 1:1 sex ratio even though the sex ratio of 
hatchling crocodilians is known to be affected by nest temperature (for reviews see 
Webb and Cooper-Preston (1989) and Deeming and Ferguson (1989)). I also assumed 
that survivorship is independent of age within a stage. This assumption does not allow 
for differential survivorship at younger ages within a stage. As more empirically 
derived data become available, a more quantitatively accurate model can be 
constructed. Matrix models have helped in understanding the relative impacts of 
management decisions on sea turtle life history stages (Crouse et al., 1987; Crowder 
et al., 1994). Crouse and co-workers (1987) and Crowder and co-workers (1994) 
showed, using elasticity analysis, that the annual survivorship of “‘large juvenile” 
loggerhead sea turtles had the greatest effect on population growth. This knowledge 
has redirected the efforts of management from focusing on headstarting to other 
critical life stages. Most crocodilians share a similar life history pattern with sea turtles 
and both have been subjected to both deliberate and incidental human impacts. It is not 
surprising to find similar matrix model results for C. acutus and those of sea turtles. 
Using the results of models like these to compare the value of alternative management 
strategies would further require an evaluation of the susceptibility of a given 
parameter to management (Nichols and Hines, 2002) and therefore the relative cost 
for a given effect (Nichols et al., 1976). In the case of crocodilians the effects of life 
history stages on population growth have not been systematically explored until now, 
although models of the type presented here are suggested by Nichols (1987). 


Simulated management scenarios—TYhe simulated management scenarios 
further strengthen the argument that sub-adult annual survivorship has the greatest 
proportional effect on population growth (see Fig. 4). These simulations demonstrate 
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Fic. 4. 20 year simulated population projections based on four management scenarios (A—C). The 
version | of the model projection is included for a baseline comparison. Simulation A increased fecundity 
50% and maintained that level of fecundity for 2 years. Simulation B increased fecundity 50% and 
maintained that level of fecundity for 20 years. Simulation C increased the annual survivorship of stage 3 
from 78% to 84% for all 20 years of the simulation. 


the relative effects on population growth of targeting particular life history stages of 
C. acutus. Increasing sub-adult annual survivorship a relatively small amount of 
2.5% (simulation C) resulted in Ag = 1.019. To approximate this level of population 
growth, a 21% increase in fecundity (simulation B) was needed. This is equivalent to 
collecting the eggs of 31.5% of the crocodile nests in Florida and incubating them in 
captivity with a survivorship to hatchling stage of 75%, as has been modeled for 
alligators (Nichols et al., 1976), before releasing them (assuming that wild and 
captive raised hatchlings have the same vital statistics). Hatch-and-release programs 
have been modeled for crocodilians and under some circumstances have been 
proposed as a viable management strategy, including for C. acutus in Florida 
(Behler, 1978; USFWS, 1984; Moler, 1991). 

Simulation A represented a more realistic hatch and release management 
scenario. This 20-year simulation considered the effects of a limited 2-year hatch 
and release program and continued to track the population for 20 years. As expected, 
simulation A tracks simulation B (Fig. 4) for the first 3 years. At year 3, the non- 
hatchling population in simulation A declined as it moved back to the stable age 
distribution of version 1, eventually increasing again at about year 8 at the same rate 
as version |. After 20 years simulation A resulted in 2.7% more non-hatchlings than 
version |. 

Moler (1991) noted that over an | 1-year period 25 C. acutus (~two per year) 
died on the two highways that connect mainland Florida to Key Largo. Similar 
mortality figures are reported by Kushlan (1988) and Kushlan and Mazzotti (1989b). 
Many of the dead C. acutus are referred to as “‘adult”’. Thus I assumed that all were 
>1 m total length and were at least sub-adults in the stage classification presented 
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here. Increasing sub-adult annual survivorship from 78.6% to 80.9% in simulation C 
was approximately equal to saving two sub-adult C. acutus in the first year of the 
simulation. Reducing road mortality may be possible by installing bridges or box 
culverts in critical areas (USFWS, 1984). Most management plans that would cause 
an increase in the annual survivorship of sub-adults would likely result in an 
increase in the annual survivorship of all stages >1 m total length. Although this 
was not explicitly modeled in simulation C, the results of the elasticity analysis 
showed that increasing annual survivorship of the adult stage would also have 
proportionally large effects on population growth (see Fig. 3). 

Although the results presented here showed that factors affecting hatchling 
survival had a small relative effect on population growth in C. acutus headstarting may 
in some cases be the only management choice. It can be of critical importance when 
trying to boost very small populations (Heppell, 1998). Headstarting can be 
problematic; for example, at Largo Enriquillo, Dominican Republic, a headstarting 
program for C. acutus failed due to personnel problems (Schubert et al., 1996). This 
population at one time was quite large (Thorbjarnarson, 1988) but has recently 
declined dramatically (Jackson, 1993; Schubert et al., 1996). It does appear that with 
increased protection, primarily by stopping illegal killing, this population may recover 
(Schubert et al., 1996), which is consistent with the results of the model presented here. 

In summary, I believe that the future of C. acutus in Florida looks hopeful. 
Version |, the best estimate of the parameters for the model, which is based to a large 
extent on published data, is consistent with observations that the Florida population 
of C. acutus is increasing. Complacency towards this population is not wise and 
hopefully the results of this modeling effort can help prioritize management efforts. 
Although they are difficult to obtain, estimates of the annual survivorship of large 
crocodilians in relation to management actions are needed. This lack of data limits 
our abilities to accurately predict the outcome of management decisions. In addition, 
the increasing human population of Florida will place increased pressure on C. acutus 
populations. Human-mediated disturbance could include road construction, increased 
road traffic and/or boat traffic, habitat destruction and Everglades restoration. 
Managers will be faced with difficult decisions and it is hoped that recognizing the 
relative effects of various life history stages, especially that of sub-adults and adults 
of C. acutus, will assist in future decision making. 
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AsstTrRact: The first documented report of escolar chino off the U.S. Atlantic coast is presented. This 
record represents a modest extension to the known geographic distribution of the species; but it is 
nevertheless a new addition to the fishes known to occur in U.S. coastal waters. Morphometric and 
meristic data are given to enhance the brief original description. The escolar chino caught off Florida had 
slightly barbed canine teeth and more scales in the lateral line than the specimen described by Poey. 
Histological analysis revealed a female with regressed (or class 1) ovaries. The taxonomy of the escolar 
chino is briefly reviewed. 


Key Words: Escolar chino, range extension, Florida 


THE escolar chino, Scombrops oculatus (Poey, 1860), (Perciformes, scom- 
bropidae) occurs in the western Atlantic Ocean on outer continental shelves and the 
slopes of seamounts (Smith-Vaniz et al., 1999). It has been characterized as being 
uncommon, demersal, bathyal, and solitary (social habits) (Claro, 1994). Escolar 
chino have been reported from Bermuda, the Bahamian side of the Straits of Florida, 
Cuba, Dominican Republic, Puerto Rico, and the U.S. Virgin Islands at depths of 
200-640 m (Poey, 1860; Zaneveld, 1983; Robins and Ray, 1986; Mochizuki, 1989; 
Claro, 1994; Silva, 1994; Smith-Vaniz et al., 1999; I.G.F.A., 2002). Robins and Ray 
(1986) suggested that this species may be distributed worldwide in tropical regions, 
but the records in the main U.S. museums indicate that this species is confined to the 
localities previously mentioned. 

We report here the first record of the escolar chino off the U.S. Atlantic coast. 
This record is based on two specimens captured 43 nautical miles (77.4 km) due east 
of Jupiter, Palm Beach County, Florida. Capt. Paul Johnson caught the fishes on the 
night of July 23, 1999, in the vicinity of 26.9400° N, 79.5900° W’ in 363-424 m of 
water. The fish were caught on hook-and-line gear rigged for commercial fishing. 
This note is based on one of the specimens because the second was cooked and eaten 
by Capt. Johnson and his crew. 


MetHops—The identification was based on the original description (Poey, 1860) and on Robins and 
Ray (1986). Although some systematists place the escolar chino in the family Pomatomidae (e.g., Robins 
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Fic. 1. Scombrops oculatus, FSBC 19124, 730 mm SL and 920 mm TL. (A) lateral view. (B) 
a typical barbed tooth (see Discussion). (C) histological section of gonad, 100, showing a gonad in 
regressed maturation stage or class 1. 


and Ray, 1986), we followed Johnson (1986) and Heemstra (1996) in assigning the genus Scombrops 
to the family Scombropidae. Meristic and morphometric data were obtained according to standard 
procedures (Hubbs and Lagler, 1949); measurements are expressed in mm and/or percentage of standard 
length (SL). A portion of the gonad was fixed in 15% buffered Formalin, embedded in plastic, sectioned to 
4 «tm, and stained with hematoxylin and eosin for histological examination. The maturation stage was 
determined according to criteria established by Taylor and co-workers (1998). The specimen was 
catalogued in the Florida Fish and Wildlife Conservation Commission Florida Marine Research Institute’s 
collection under catalogue number FSBC 19124. 


ResuLts—Scombrops oculatus, FSBC 19124.—Diagnostic characters: Body 
elongate, head slender; soft dorsal and anal similar in shape and height, their outlines 
concave (Fig. 1A); first dorsal with 8 spines, anal with 3 spines; mouth very large, 
lower jaw projecting, maxillary extending beyond pupil. Teeth in each jaw in a sin- 
gle row; canines slightly barbed, long, sharp, well separated; upper jaw with 2 
prominent fang-like teeth anterior to the vomer and with 11 large canines on 
premaxillary; lower jaw with 2 small canines at front, followed by 9 large canines. 
Teeth on vomer and palatines in a single row; vomerine, lingual, and pharyngeal 
teeth villiform. Preopercle entire; opercle with 2 weak spines. Pectoral and pelvic 
fins short; pelvics inserted under the pectorals. Caudal fin forked. Eyes large. Scales 
large, cycloid, deciduous, with entire edges. 
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Meristic characters: D VIII, I-14; A H,12; P; 15; P> 1,5. Scales in lateral line 60. 
Gill rakers on Ist arch: 17 (1 + 5 rudiments upper; 7 + 4 rudiments lower). The 
underside of both opercula has a pseudobranch with 73 free lamellae. There are 5 
complete gill arches (holobranchs). 

Morphometric characters (in mm and % SL unless otherwise indicated): SL 
730; total length (TL) 920; body depth 210 (28.7%); head length (HL) 226.5 (31%). 
Orbit diameter 64.1 (28.3% HL); interorbital distance (width) 69.5. Predorsal 
distance 266 (36.4%). Spinous dorsal fin base 130.4 (17.8%); pectoral fin length 
133.9; pectoral fin base 36.6; pelvic fin length 107.7; pelvic fin base 93.6; caudal fin 
length 178.0; caudal peduncle depth 69.0. Longest dorsal fin spine (2nd) 105.2; anal 
fin height 106. Upper jaw length 133.4; lower jaw length 113.3. 

The weight of the preserved specimen is 16.82 kg (37 Ibs). 

Color: The natural coloration of the specimen after 24 hrs on ice was violet 
brown, paler ventrally. The iris was golden orange. The preserved fish is dark brown. 

Gonad histology: The analysis revealed that this female fish was in regressed 
maturation stage or Class | on July 23 (Fig. IC). 


DiscussioN—The escolar chino was originally described as Latebrus oculatus by 
Poey in 1860. Giinther (1887) assigned Latebrus to subgeneric rank under the genus 
Scombrops. Ginther (1887) remarked that, ““but for the barbed teeth, he would be 
disposed to unite S. (Latebrus) oculatus with S. chilodipteroides,” a species found in 
the Pacific Ocean. Goode and Bean (1896) pointed out that they would not follow 
Giinther (1887) because the scales of S. chilodipteroides are considerably smaller, 
the body is more slender, the fins are longer, the two dorsals are more closely 
approximated, the lower jaw appears to be longer, and the second dorsal is placed 
much farther forward, relative to the vent. Jordan and Evermann (1896) also noticed 
that in Poey’s figure, the premaxillary canines were represented as distinctly barbed, 
although nothing was said of this character in the description. Most importantly, 
Jordan and Evermann (1896) suggested that the presence of barbed teeth may per- 
haps define the subgenus or genus Latebrus, because the teeth in the Japanese 
S. chilodipteroides are not barbed. 

The specimen of escolar chino collected in Florida has slightly barbed canine 
teeth (Fig. 1B); this character is clear and magnified by color in the fresh specimen. 
The slightly barbed (distal) portion of the teeth was dark yellow and the lower 
(proximal) portion white. There are more scales in the lateral line of FSBC 19124 
than in the lateral line of the specimen described by Poey (1860). A possible cause 
for this difference is that the scales in the transverse series are rather irregularly 
arranged and do not correspond with the number of scales along the lateral line. A 
report on the taxonomic implications of these findings is in preparation. 

The escolar chino is also known by the common names Atlantic scombrops 
(Robins and Ray, 1986) and boral prieto in the Dominican Republic (Silva, 1994). 
The common name pseudosnapper will appear in the upcoming FAO guide to 
western Atlantic fishes (Gilbert, 2003). The escolar chino caught off Florida rep- 
resents a modest extension to the known geographic distribution of the species. It is 
nevertheless a new addition to the fishes known to occur in U.S. coastal waters. 
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ABSTRACT: Of the 18 species of sunflowers recorded from Florida, Helianthus carnosus Small is 
easily distinguished by the scapose flower stalks, yellow disk flowers, and glabrous leaves. This perennial 
herb is endemic to northeastern Florida and is considered endangered by the Florida Department of 
Agriculture and Consumer Services. The lakeside sunflower is reported for the first time from Putnam 
County representing the only known population west of the St. Johns River. Associated plants are listed, 
and the habitat is described. 


Key Words: Conservation, rare plant, listed species, flatwoods 


On August 20, 2001 we found three individuals of Helianthus carnosus 
blooming at the edge of a mown trail at the St. Johns River Water Management 
District Natural Area west of the Kay Larkin Airport in Palatka. The coordinates of 
the site were measured with a hand-held Garmin III Plus GPS receiver as north 
latitude 29° 40.047’ and west longitude 81° 41.877’. The SJRWMD Natural Area is 
an 8.8-acre parcel consisting mostly of flatwoods. The forest canopy was a mixture 
of mature longleaf pine (Pinus palustris) and slash pine (Pinus elliottii) with an 
understory dominated by saw palmetto (Serenoa repens) and gallberry (Ulex glabra). 
The Helianthus carnosus plants were located near a grassy depression containing a 
few clusters of hooded pitcherplant (Sarracenia minor). This site was intermittently 
grazed until 1991, when the District purchased the parcel, but has not burned 
for many years. 

A voucher specimen of Helianthus carnosus was collected and deposited in the 
herbarium at the University of South Florida, Institute for Systematic Botany, 
Tampa, Florida (USA). The two remaining plants flowered through the end of 
August. By 25 October 2001, the distal portions of the flower stalks were dead and 
dried. Chafin (2000) listed flowering as occurring from August through October. 
Most of the year, lakeside sunflower exists as a rosette of leaves at the ground 
surface. From our observations, Helianthus carnosus tolerates mowing well, which 
may help the plant survive at sites where fire has long been suppressed. 

The habitat of lakeside sunflower is wet prairies and flatwoods (Wunderlin, 
1998; Chafin, 2000; Coile, 2000). We investigated the vegetation associated with the 
Helianthus carnosus plants on 29 October 2001. In the vicinity of the three original 
plants that had flowered, we located three, four, and seven rosettes of Helianthus 
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TABLE 1. Plants found in one square meter quadrats centered on three different clusters of 
Helianthus carnosus rosettes from Putnam Co., Florida. An ‘X’ indicates presence within the quadrat. 


Quadrat 
Family Species 1 2 3 
Monocotyledons 
Hyopoxidaceae Hypoxis juncea Sm. x — — 
Poaceae Andropogon brachystachyus Chapm. — x — 
Poaceae Aristida sp. xX — — 
Poaceae Dichanthelium ensifolium (Baldwin ex Elliott) Gould x x xX 
Poaceae Dichanthelium strigosum (Muhl. ex Elliott) Freckmann x x xX 
Poaceae Paspalum setaceum Michx. x — x 
Xyridaceae Xyris ambigua (Beyr. ex Kunth) x — — 
Dicotyledons 
Aquifoliaceae Ilex glabra (L.) A. Gray — x XxX 
Asteraceae Elephantopus elatus Bertol. x x XxX 
Asteraceae Pityopsis graminifolia (Michx.) Nutt. — x — 
Campanulaceae Lobelia paludosa Nutt. x x o 
Cistaceae Lechea minor L. x _— — 
Clusiaceae Hypericum brachyphyllum (Spach) Steud. — XxX — 
Ericaceae Gaylussacia dumosa (Andrews) Torr. & A. Gray — Xx x 
Ericaceae Lyonia fruticosa (Michx.) G. S. Torr. — x — 
Ericaceae Vaccinium myrsinites Lam. — — xX 
Fabaceae Mimosa quadrivalvis L. x — — 
Fagaceae Quercus minima (Sarg.) Small xX » xX 
Loganiaceae Gelsemium sempervirens (L.) W. T. Aiton — xX — 
Rubiaceae Hedyotis procumbens (J. F. Gmel.) Fosberg — x — 
Rubiaceae Hedyotis uniflora (L.) Lam. — x — 
Scrophulariaceae Gratiola hispida (Benth. ex. Lindl.) Pollard xX — xX 
Totals WW 14 2) 


carnosus. These rosettes appeared to be derived from rhizomes rather than seeds. 
Within one square meter quadrats centered on the three patches of Helianthus 
carnosus rosettes, we found 22 species of vascular plants representing 14 families 
(Table 1). Grasses had the highest cover, visually estimated to range between 50— 
80% within the quadrats. The vegetation adjacent to the mown trail was domi- 
nated by dense shrubs, especially gallberry and saw palmetto. Other nearby plants 
included longleaf pine, slash pine, water oak (Quercus nigra), fourpetal st. john’s- 
wort (Hypericum tetrapetalum), muscadine grape (Vitis rotundifolia), and earleaf 
greenbrier (Smilax auriculata). 

The soil supporting the Helianthus carnosus population was Pomona fine sand, 
a sandy, siliceous, hyperthermic ultic haplaquod (Baldwin, 2001). This spodosol was 
poorly drained and had a sandy clay loam subsoil (Readle, 1990). Water levels in 
a nearby shallow well fluctuated from ground surface to 1.3 m below the surface in 
weekly measurements taken from 1995 to 1999. Average seasonal high water levels 
in the well occurred in January, February, and March. The lowest levels were in May, 
June, and July at the end of the dry season into the beginning of the rainy season in 
northern Florida. 
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Prior to this discovery, only nine populations of Helianthus carnosus were 
known in St. Johns, Flagler, and Volusia counties, Florida (Chafin, 2000). Although 
this rare plant was previously recorded from Clay County (Wunderlin and Hansen, 
2002), it has not been observed west of the St. Johns River for many decades (Chafin, 
2000). This Putnam County site is significant because it is the only known popula- 
tion on publicly owned land. Furthermore, it raises the possibility that Helianthus 
carnosus may be more widely distributed in northeastern Florida, but has been 
overlooked. 
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ABSTRACT: Infection with upper respiratory tract disease (URTD) has been linked to population 
declines of tortoises. Here we test whether gopher tortoises (Gopherus polyphemus) at Archbold 
Biological Station (ABS), a long-protected area in south-central Florida, have been exposed to, or are 
currently infected with Mycoplasma agassizii, a causal agent of URTD. We examined all gopher tortoises 
encountered for clinical signs of URTD. In addition, we collected blood (n= 11) and nasal samples 
(n = 10) from tortoises found throughout ABS. We used PCR-based tests of nasal samples to determine if 
gopher tortoises are currently infected with M. agassizii and ELISAs for blood samples to test if exposure 
to this mycoplasmid has taken place. No tortoise at ABS showed any clinical signs of URTD. 
Furthermore, all blood samples were seronegative, indicating that no tortoise had developed a detectable 
immune response to M. agassizii. All nasal lavages were also negative for the presence of M. agassizii. 
Thus, gopher tortoises at ABS have not been exposed to this mycoplasmid. Comparisons of characteristics 
between populations with varying levels of URTD prevalence, including sites such as ABS without any 
cases of infection, will help decipher the role that this disease plays in population declines of gopher 
tortoises. 


Key Words: Gopherus, Upper Respiratory Tract Disease, Mycoplasma, Archbold 
Biological Station 


GOPHER tortoises (Gopherus polyphemus) are significant contributors to their 
xeric ecosystems within the southeastern United States. The most fossorial of the four 
tortoises found in North America, gopher tortoises dig burrows averaging 4.5 m long 
and 2 m deep (Hansen, 1953; Diemer, 1986; Kaczor and Hartnett, 1989; Brown et al., 
1999). Nearly 400 species of invertebrates and vertebrates have been observed using 
gopher tortoise burrows. In fact, many have recognized gopher tortoises as a keystone 
species (Brown et al., 1999). Thus, the health of gopher tortoise populations may 
have a large impact on other species in the southeastern United States. 

Recently, an upper respiratory tract disease (URTD) has been found in 
populations of tortoises. Though originally discovered in desert tortoises (Xerobates 
agassizii), the disease has been found in other species as well (Jacobson, 1991). One 
of the causal agents, Mycoplasma agassizii, is a virulent bacterium that resides in 
the respiratory tract of infected tortoises (Brown et al., 1994, 2001). This disease is 
characterized by mild to severe nasal and ocular discharge—though a tortoise can be 
infected without showing clinical signs, and may remain a carrier for life (Jacobson, 
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1993). Upper respiratory tract disease is horizontally transmitted through direct 
contact, although airborne transmission is possible over short distances (<0.5 m; 
Jacobson et al., 1995). 

Presence of URTD has been linked to population declines in tortoises (Jacobson, 
1993: Corn, 1994; Brown et al., 2001; Gates et al., 2002). Whether the disease itself 
causes such declines is unclear, however infection with URTD 1s capable of killing 
tortoises (Jacobson, 1993). Understanding the link between URTD and popula- 
tion declines requires 1) identifying populations with and without URTD and 2) 
monitoring populations that differ in URTD prevalence over time. Another reason to 
evaluate the distribution of URTD is to limit the possibility of human-induced 
spreading of the disease between populations. Gopher tortoises are a protected spe- 
cies. Accordingly, gopher tortoise relocations are often necessary for conservation 
purposes. However, translocations of gopher tortoises without prior testing for URTD 
may have aided in the spread of M. agassizii among populations (Diemer, 1987). 
Upper respiratory tract disease has no cure; therefore, determining the distribution of 
the disease is necessary in order to make informed decisions regarding conservation 
initiatives and proper land management practices. 

The distribution and prevalence of URTD in gopher tortoises were recently 
tested at 53 sites in Florida (Diemer Berish et al., 2000). However, Diemer Berish 
and coworkers (2000) did not sample the gopher tortoise population at Archbold 
Biological Station (ABS), Highlands County, Florida. ABS is a 5,140-acre preserve 
located at the southern end of the Lake Wales Ridge (Swain, 1998). The station 
houses a large population of gopher tortoises, with ~ 1200 individuals marked since 
1967. Many tortoises have been relocated to ABS, both with and without permission 
(Layne, 1989). Known relocations have provided an opportunity to study the effects 
of relocation on behavior and life history (Layne, 1989). However, such relocations 
may have inadvertently introduced pathogens such as M. agassizii to gopher tortoises 
at ABS—if they were not already present. Here we test for presence of M. agassizii 
at ABS. The status of infection at ABS is critical to future management of this 
population. Further, once the prevalence of URTD at ABS is known, the population 
can be used in comparative studies to evaluate the relationship between URTD and 
population viability. 


MetHops—Presence of clinical signs, such as mucous discharge from the nasal and/or ocular region, 
in conjunction with positive ELISA and PCR test results, indicates infection with URTD (Brown et al., 
1999; Jacobson et al., 1995; Diemer Berish et al., 2000). An ELISA test uses blood samples to detect 
antibodies to M. agassizii, thus indicating that a tortoise has been exposed to the bacterium at some point in 
its life (Schumacher et al., 1993). Evaluating nasal lavage samples using a PCR-based test determines 
whether a tortoise is currently infected with M. agassizii (Brown et al., 1995). Should positive results arise 
from both the ELISA and PCR analysis without clinical signs, it is possible that the tortoise has subclinical 
infection or is colonized with a nonpathogenic strain. Currently no tools are available to designate a strain 
as pathogenic or not, thus we employ all three methods to detect potential cases of URTD. 

From September 2001—July 2002 we surveyed ABS for gopher tortoises. Each individual 
encountered was measured and marked following Layne (1989), and examined for external signs of 
URTD. Such signs include 1) nasal discharge, 2) ocular discharge, 3) conjunctivitis, and 4) swollen 
eyelids (Brown et al., 1999; Jacobson, 1993). 

We transferred 14 gopher tortoises to a laboratory to collect blood and nasal samples. Blood was 
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TABLE 1. Evidence of URTD for gopher tortoises at Archbold Biological Station. Results are 
presented as the percentage of gopher tortoises testing positive or negative for clinical signs of URTD or 
infection with M. agassizii (ELISA and PCR). N = number of individuals tested. 


Test N % Positive % Negative 
External symptoms 103 0 100 
ELISA 11 0 100 
PCR 10 0 100 


drawn (0.3—1.0 ml) from the brachial vein of 11 tortoises and placed into tubes containing lithium heparin. 
We then centrifuged all samples and pipetted the plasma into transport tubes. In addition, we flushed each 
naris of 10 tortoises with 0.5 ml of sterile phosphate-buffered saline to collect nasal lavage samples. All 
samples were frozen and transported to the Mycoplasma Testing Lab at the University of Florida, 
Gainesville, FL for diagnostic tests (ELISA and PCR-based). 


RESULTS—We examined 103 gopher tortoises for clinical signs of URTD, 
including 46 adult females and 37 adult males (20 tortoises lacked any sexually 
distinguishable characteristics). No tortoise showed any clinical signs of URTD. 
That is, no tortoise had mucous discharge from the nose or eyes, conjunctivitis or 
swollen eyelids (Table 1). Some individuals had a small amount of clear discharge 
from the corners of the eyes and/or naris, but it was temporary and onset, duration, 
and completion were associated with periods of manual handling. 

We tested 11 individuals (7 females, 3 males, 1 unidentified sex) for presence of 
antibodies to M. agassizii (ELISA test). These tortoises ranged from 1.9—7.8 kg 
(mean = 4.1 + 1.8 kg SD) mass and 232—343 mm (mean = 279 + 35 mm SD) 
straight-line carapace length. One hundred percent of the samples were negative for 
antibodies to M. agassizii (Table 1). 

We tested ten tortoises (4 females, 5 males, 1 unidentified sex) for current 
infection by M. agassizii using nasal samples and the PCR-based test. These 
tortoises ranged from 1.9-5.7 kg (mean = 3.5 + 1 kg SD) mass and 232-312 mm 
(mean = 266 + 23 mm SD) straight-line carapace length. Eight of these tortoises 
were also used for ELISA testing. All samples were negative for presence of 
M. agassizii (Table 1). 


DiscussioN—No gopher tortoise at ABS showed any clinical signs of URTD. 
Further, no individual tested positive for exposure to M. agassizii using the ELISA 
or PCR test (Table 1). Because we sampled gopher tortoises from throughout ABS 
property, we conclude that URTD is most likely not present at ABS. However, the 
small sample size used in our study leaves open the possibility that URTD is present 
but went undetected. Further testing is needed to definitively rule out presence of 
URTD. 

Our results increase the list of populations in Florida tested for URTD. Diemer 
Berish and coworkers (2000) found 74% of 53 sites were negative for URTD. How- 
ever, the site surveyed closest to ABS, Platt Branch Mitigation Park (Highlands 
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County; 16 km from ABS), had 17/20 individuals test positive for URTD. While 
gopher tortoise home ranges are typically less than a hectare (Diemer, 1992b), a 
strong possibility of the disease being introduced to ABS exists due to the high 
occurrence of local gopher tortoise translocations. These translocations are mainly 
the result of a tremendous population growth in Highlands County, with a 28% 
increase in population size from 1990 to 2000, as well as a 163% increase in private, 
non-farm establishments from 1990-1999 (U.S. Census Bureau, 2003). Protected 
status of gopher tortoises in Florida often necessitates such translocations as 
a conservation alternative because such a tremendous growth in economic develop- 
ment is associated with a decrease in suitable gopher tortoise habitat (Diemer, 1987; 
1992a). 

As long as translocations remain the primary conservation procedure for 
protecting gopher tortoises in Florida, the status of URTD within the population 
being relocated and the receiving population must be taken into account (Diemer, 
1987). While the exact nature of how M. agassizii becomes active within a popu- 
lation is still unknown, recent evidence suggests it is due to an infected tortoise 
being introduced into a previously uninfected population, followed by horizontal 
transmission between adult tortoises (Jacobson et al., 1995; Smith et al., 1998)— 
though predisposing factors have also been implicated in abetting the development 
and spread of the disease in free-ranging populations (Jacobson et al., 1991). Thus, it 
is possible that 1) an infected tortoise has yet to be introduced to ABS despite the 
numerous translocations that have taken place or 2) the tortoises at ABS maintain 
a predisposition that prevents the development and/or the spreading of the disease. 

Further knowledge of URTD and research on its effects within tortoise 
populations is necessary in order to apply appropriate conservation tactics and 
intervention procedures. Obviously, preventative measures will sustain negative 
URTD status for a population, and until a thorough understanding of URTD is 
achieved, all efforts should be put into maintaining such practices (1.e., extensive 
testing prior to and after relocations of both participating populations; follow-up 
testing for maintenance purposes). Comparing tortoises at ABS with those in areas 
infected with M. agassizii, for instance at Platt Branch Mitigation Park, can help 
clarify how presence of the pathogen influences populations of gopher tortoises. 
Particular attention should focus on whether the spread of M. agassizii and, 
ultimately, URTD, is related to environmental factors (e.g., population density) or 
genetic factors (e.g., some individuals potentially immune to infection by 
M. agassizii). 
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AsstTRraAct: Surveillance for gastrointestinal parasites among Mayan Indians was conducted in San 
Antonio Palopo, Santa Catarina Palopo, Santa Cruz Laguna and Jaibalito, four villages around Lake 
Atitlan, Guatemala. In January and June 2001, 234 patients with gastrointestinal symptoms of abdominal 
pain were identified and their stools examined for parasites. Feces were collected from patients and zinc 
sulfate flotation was performed on each specimen. Both helminth and protozoan parasites were detected 
in varying rates from people in the different villages. Ascaris lumbricoides was the most frequent helminth 
infection detected followed by Trichuris trichiura and Hymenolepis nana. Entamoeba histolytica/dispar 
was the most common protozoan infection followed by Giardia duodenalis. Results demonstrated that 
parasite transmission is largely in children and the prevalence rate of infection is related to the time of 
year. This identification of the prevalence of parasites in different villages in relation to the current risk 
factors is suggestive of the type of preventative measures necessary to reduce parasitic infection rates. 


Key Words: Ascaris, Entamoeba, Trichuris, Hymenolepis, Guatemala 


HUMAN parasitic diseases are of public health importance accounting for 10% of 
the cases of gastrointestinal diseases in developing countries (Crompton, 1999). Of 
these, some of the highest infection rates in the world are reported in Guatemala. 
Ascaris lumbricoides has been found to have a prevalence of 15-74% in adults 
and 28-89% in children (Anderson et al., 1993; Farthing et al., 1986; Grazioso et al., 
1993; Gupta and Urrutia, 1982; Melvin and Mata, 1971; Villar et al., 1989). 
In addition, prevalence rates for Giardia duodenalis are 2—14% and 15-22%, 
Entamoeba histolytica/dispar, 3—33% and 4-39%, and Trichuris trichiura 4-64% 
and 25-69%, respectively for adults and children. A 1.0% prevalence rate for 
Hymenolepis nana in Guatemalan adults is reported. 

Recent publications have not assessed the prevalence of the various parasitic 
infections among symptomatic people living in the culturally isolated villages 
around Lake Atitlan. Therefore, our laboratory surveillance of gastrointestinal para- 
sites was completed in order to identify symptomatic cases and provide data that 
could link parasitic diseases to specific risk factors associated with living in this 
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geographical region of Guatemala. The study was conducted in January and June 
2001 to survey the prevalence of human parasites in areas with high risk factors. The 
purpose was to determine if the age or sex of the patients or if the time of year 
correlated with differences in the intensities of symptomatic parasitic infections in 
people living around Lake Atitlan. 


PaTIENTS/METHODS—Area and population—Several villages around Lake Atitlan were chosen to do 
the study. Each location was suspected of having high rates of parasitic infections and a population that 
was representative of the geographic area. Of the selected communities, there were 9,024 residents in San 
Antonio Palopo, 959 residents in Santa Catarina Palopo, 3,206 residents in Santa Cruz La Laguna and 376 
residents in Jaibalito in 1995 (De Paz Jaenike, 1997). The villages were located 40 km south of Guatemala 
City in a volcanic mountainous area. The climate was equatorial with a long rainy season that lasts from 
April to mid November. Most of the residents were engaged in subsistence, in the making of arts/crafts to 
sell, or in the farming of beans, corn, onions, squash and coffee, all on the sides of the mountains. 


Subjects—A total of 234 symptomatic patients, ages 1 to 75 years (56% females, 44% males, 75% 
children and 25% adults) from four villages around Lake Atitlan were examined during visits on 15—20 
January and 11—16 June in 2001. The median patient age was 15.3 years (range, 1-73 years), with 9.8 
years (range, 1-73) for males and 22.0 years (range, 1-67) for females. The majority of the patients 
examined were Mayan Indians. 

The objectives of the study were explained to members of the community, to their leaders and to 
school and health administrators. Symptomatic patients typically presented with diarrhea and abdominal 
symptoms, including cramps, loss of appetite, pain, bloating or tenderness. All infected adult patients were 
treated and each infected child aged | to 18 years was treated while in the presence of one of his or her 
parents in our temporary clinic area. The study was reviewed and approved by the Institutional Review 
Board at Nova Southeastern University. Informed consent was obtained from subjects, parents and/or 
guardians. 


Specimen collection—NSU-DOCARE medical teams at each site collected fecal samples directly 
from the anus of each patient. The method of collection was performed by inserting a lubricated, gloved 
finger into the rectum of patients followed by collection of feces from the glove. Each specimen was 
assigned a number that correlated with that individual patient. Investigators working on parasitic infections 
in endemic areas have not previously reported on the effectiveness and practicability of this approach. 


Laboratory analysis—¥ecal material was washed from the glove using a zinc sulfate solution and the 
fecal solution was allowed to set in a 10.0-ml glass vial for 15 minutes allowing the parasitic life stages to 
float to the top and be collected on a plastic cover slip. Cover slips were mounted onto glass slides along 
with a drop of Lugol’s iodine and observed at 100 and 400X magnification using a bright field microscope. 
All identified parasitic infections were recorded for each patient. 


Data collection—The age and gender of each patient was recorded at the time they were examined. 
Differences between group proportions (positive or negative for a parasite) were compared using the xe 
test or McNemar’s test. P values less than 0.05 were considered to be significant. 


ResuLts—Prevalence of parasites in the four villages—Of the total patients 
from all villages, the overall prevalence of intestinal parasite infections ranged 
between 25% and 47%, with a median of 38%. The most common parasite found was 
A. lumbricoides and the least common was G. duodenalis. Patients positive for A. 
lumbricoides were also positive for T. trichiura on three occasions and with H. nana, 
five times. This is likely the result of a free interaction between the contaminated 
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TABLE |. Prevalence of parasitic infections of symptomatic residents* living in villages around Lake 
Atitlan, 2001. 


No. (%) Prevalence 


San Antonio Santa Catarina Santa Cruz 
Palopo Palopo Laguna Jaibalito 
Parasitic infection n = 108 n = 89 n = 24 n= 13 
Entamoeba histolytica/dispar 2 (1.9) 5 (5.6)° 0 (0) 0 (0) 
Giardia duodenalis 1 (0.9) 4 (4.5)? 0 (0) 0 (0) 
Hymenolepis nana 3 (2.8)° 4 (4.5)° 0 (0) 1 (7.7) 
Ascaris lumbricoides 28 (25.9)? 25 (28.1)° 3 (1.3) 3 (23.1) 
Trichuris trichiura 4 3.7) 4 (4.5) 37155) 2 (15.4) 
Overall rate 38 (35.2) 42 (47.2) 6 (25.0) 6 (46.2) 


* Patients with gastrointestinal complaints. 
> significantly different from Santa Catarina Palopo and Jaibalito. 


environment and susceptible people because results showed that individuals infected 
with A. lumbricoides did not have a greater likelihood of being infected with 
T. trichiura or H. nana than those who were not infected. 

San Antonio Palopo and Santa Catarina Palopo had a significantly higher 
prevalence of E. histolytica/dispar infection than Santa Cruz La Laguna or Jaibalito 
(Table 1). H. nana infections were not found in Santa Cruz La Laguna and there was 
also a marked absence of E. histolytica/dispar and G. duodenalis in Santa Cruz 
La Laguna and Jaibalito. Santa Cruz La Laguna had a lower infection rate for 
A. lumbricoides and Jaibalito had a higher infection rate for T. trichiura than the 
other villages. However, the overall rate of parasitic infections was similar for all 
four villages. 


Seasonal effect—Of the four villages examined for parasitic infections, San 
Antonio Palopo and Santa Catarina Palopo had large numbers of patients and were 
used to determine seasonal prevalence (Table 2). Fewer infections with E. histolytica/ 
dispar, G. duodenalis, H. nana and A. lumbricoides and an overall lower prevalence 
rate of parasitic infections was seen in June than in January in San Antonio Palopo. 
E. histolytica/dispar was not found in June in Santa Catarina Palopo. 


Age—The children in each of the four towns had more infections than the adults 
did, with A. lumbricoides being most prevalent in children (Table 3). Examination of 
175 children revealed an increased prevalence rate for infections with E. histolytica/ 
dispar, G. duodenalis, H. nana and A. lumbricoides in children living in San 
Antonio Palopo and Santa Catarina Palopo compared to Santa Cruz La Laguna and 
Jaibalito. Adults were found to have a lower prevalence rate of infections when 
compared to children. 


Sex—Host gender had a pronounced influence on the prevalence of infections 
(Table 4). A. lumbricoides was more common in female than in male children and E. 
histolytica/dispar was present in females but absent in male adults. Overall, more 
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TABLE 2. Seasonal prevalence of parasitic infections of symptomatic residents* of San Antonio 
Palopo and Santa Catarina Palopo, 2001. 


No. (%) Prevalence 


San Antonio Palopo Santa Catarina Palopo 
January June January June 
Parasitic infection n= 84 n= 24 in == Syl fi— 0 
Entamoeba histolytica/dispar 2 (2.0) 0 (0)? 5 (9.8) 0 (0)? 
Giardia duodenalis 1 (1.0) 0 (0)? 1 (2.0) 3 (7.9) 
Hymenolepis nana 3 (3.0) 0 (0)? 1 (2.0) 3 (7.9) 
Ascaris lumbricoides 28 (28.0) 0 (0)? 16 (31.4) 9 (23.7) 
Trichuris trichiura 3 (3.0) 1 (4.2) 1 (2.0) 3 (7.9) 
Overall rate 37 (37.0) 1 (4.2) 24 (47.1) 17 (44.7) 


* Patients with gastrointestinal complaints. 
> Significantly fewer infections in June than in January. 


infections were found in female patients (24.4%) than in male patients (15.8%) with 
A. lumbricoides being most prevalent in females. Female children had more G. 
duodenalis, H. nana and A. lumbricoides infections than female adults and male 
children had more parasitic infections than male adults did. 


DiscussiIoN—This is the first reported study demonstrating prevalence rates 
for gastrointestinal parasitic infections in villages around Lake Atitlan. The high 
prevalence of parasitic infections in these four culturally isolated villages was lower 
than previous reports from distant towns using different collection techniques 
(Anderson et al., 1993; Farthing et al., 1986; Gupta and Urrutia, 1982; Mata, 1982; 
Melvin and Mata, 1971). The observed decline of parasitic infections might have 
been the result of attempts by others to treat patients. 

Parasitic infections were associated with several specific risk factors. First, the 
Mayans in the rural areas did not obtain medical care on a regular basis due to 
a shortage of healthcare workers and the expense of medical care. Second, the 
villagers were exposed to infection by using Lake Atitlan for drinking water, laundry 
of clothes, bathing and fishing because the lake did not have a river feeding it or 
leaving it, therefore collecting everything that enters it. Third, there was consider- 
able run-off of rainwater into the lake and the water contained organic matter from 
the crops, livestock and humans living in the villages. 

The prevalence of parasitic infections and transmission rates were compared 
between the different but closely related communities. Transmission was not similar 
year around for E. histolytica/dispar, G. duodenalis, H. nana, and A. lumbricoides, 
but was for T. trichiura, a distinction that is most evident during times when 
transmission is at a relatively low level. Differences in the prevalence rates between 
villages could not be explained by differences in population densities, environmental 
features or personal habits of the residents. 

Seasonal differences in prevalence rates coincided with the periods of extremely 
dry climate (December to April) and drenching wet months (May to November). 
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Fewer infections in June followed the dry season and higher infections in January 
followed the wet season. Possibly, dry weather inhibited the survival of parasite life 
forms, killing them by desiccation in the soil and the wet season supported the 
survival and transmission of the parasites. 

Our data demonstrate that parasitic infections around Lake Atitlan occur early 
in life and are more frequent in young males and females than in adults. We found 
that children, age 1 to 18, suffered the highest rates of infection and the most ill- 
ness. Possibly this is because children are the most sensitive indicators of parasite 
transmission. This decrease in prevalence of infection with age was consistent with 
the findings that there was a noticeable decrease in infection with parasites in people 
over 28 years of age (Anderson et al., 1993; Farthing et al., 1986). An obvious 
means by which the decreasing rate of infection with age may be caused is through 
a gradual decrease in the proportion of the aging population that is susceptible to 
infection. In addition, many of the men in the communities may have lower infection 
rates because they spend many hours working outside the downtown, contaminated 
areas each day. The areas around the schools are contaminated with human feces and 
the children are regularly exposed to parasites. 

We found differences in parasitic infection rates due to the patient’s age and/or 
sex and found that infection rates were related to the time of year the study was 
done. Evidence implies that there are significant numbers of infected individuals 
with or without symptoms that could lead to larger outbreaks of parasites in the local 
villages. Such infections pose a substantial health threat in these communities and 
measures to improve environmental sanitation in the villages and lake are needed. In 
conjunction with this, education of the residents on methods of maintaining adequate 
personal hygiene would result in a coordinated approach to prevention and control 
of parasitic infections around Lake Atitlan. 
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ABSTRACT: Studies were conducted to identify the behavioral components of maternal behavior of 
the shrew, Blarina carolinensis, in non-manipulated litters, and to assess the effects of litter gender 
composition (LGC) on subsequent maternal behavior. Maternal behaviors consisted of nursing, licking, 
nest maintenace, retrieval of young, and miscellaneous (non-maternal) activities including eating, 
drinking, grooming, scratching, and sleeping. Weaning occured between 23-24 days of age. On day 13, 
when male and female pups could first be distinguished from one another on the basis of the appearance 
of mammae in females, 13 litters containing 5 pups each were culled to 4 and randomly assigned to I of 3 
treatment groups: MM (all males), MF (2 males and 2 females), and FF (all females). Behavioral 
interactions between mothers and their young were observed from day 2 to day 22. There were no 
significant differences in mass between the treatment groups for individual neonates and litters over the 
entire weaning period. At birth, neonates ranged from 0.67-0.74 g, and the weight of the entire litters 
ranged from 2.43—3.17 g. At weaning, the pups weighed between §.84—9.47 g, with no significant gender 
differences. Maternal behavior was more pronounced in the days immediately following birth. There was 
a significant decrease in the frequency of nursing and licking between days 2 and 22 following birth, and 
an increase in miscellaneous activities. The amount of time spent in nursing, licking, and nest maintenance 
decreased significantly, with a similar increase for miscellaneous activities. MM mothers showed a higher 
frequency and duration of nursing and licking than those from the MF and FF groups. In B. carolinensis, 
LGC affects some aspects of maternal behavior but not the physical development of offspring. 


Key Words: Blarina carolinensis, litter gender composition, maternal behavior 


PREVIOUS research on mice (Mus musculus) has shown that litter gender 
composition (LGC) can affect various aspects of maternal behavior directed toward 
pups (Alleva et al., 1989), as well as the subsequent exploratory behavior and emo- 
tionality (Denenberg and Morton, 1964), agonistic interactions (Mendl and Paul, 
1991), reproduction (Birke and Sadler, 1987), and central nervous system develop- 
ment (Laviola and Alleva, 1995) of the pups themselves. Similar findings have been 
reported for Norway rats, Rattus norvegicus (Sharpe, 1975; De Cabo et al., 1992), 
house cats (Caro, 1981), and several species of primates (Barrett et al., 1992). The 
possible effects of LGC on other mammalian taxa have received little attention, 
including shrews (Churchfield, 1990; Laviola and Alleva, 1995). In addition, there is 
little detailed information on maternal behavior in shrews (Churchfield, 1990). 

The southern short-tailed shrew (Blarina carolinensis) ranges from southern 
Nebraska eastward throughout most of the southeastern U.S. (Whitaker, 1996). In 
Florida, it typically inhabits small burrows under plant debris and decaying logs 
(Genoways and Choate, 1998), and prefers mesic habitats in bayheads and flatwoods 
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(Gingerich, 1999). The only detailed study on reproduction in B. carolinensis was 
conducted on a population from Aiken, SC, and reported litter sizes ranging from 2— 
6, with a mean of 3.75 (O’Farrell et al., 1977). Previous success at captive breeding, 
as well as analyses of the interactions between B. carolinensis and several species of 
sympatric rodents (Punzo, 2003), gave me the opportunity to study other behavioral 
parameters in this shrew. The purpose of the present study was to identify the 
behavioral components of maternal behavior in B. carolinensis and the effects of 
LGC on the interactions between the maternal parent and her offspring. To my 
knowledge, no data exist on potential consequences of LGC on behavior for any 
species of shrew. 


METHOoODs—AIl of the shrews used in these experiments were obtained from a captive-bred colony 
that was initiated from adult shrews originally collected in Leon Co., FL, in 1998. Adult shrews ranged in 
weight from 14.3—21.2 g (mean: 18.1 g). The procedures used to house shrews have been described 
elsewhere (Punzo, 2003). To summarize, 10 pairs of male and female shrews were housed in 30 X 18 x 15 
cm plastic rodent cages (one pair/cage) and provided ad libitum with food and water. Their diet consisted 
of an equal mixture (by weight) of commercial cat chow (Ralston Purina, St. Louis, MO) and earthworms 
(Lumbricus terrestris). The shrews were provided with 15 g of freshly prepared food daily. Animals were 
maintained at 22° + 1.0°C, 60-65% RH, and a 10:14 white light/red light cycle (red lights [dark phase] 
turned on at 0900 hr). 

Ten randomly selected pregnant females, identified through vaginal smears, were placed individually 
in 1.0 X 1.0 X 0.5 m transparent plastic arenas. The bottom of each arena consisted of clear plexiglass so 
that observations could be made from above and below the arena. The arena floors were provided with a 5.0- 
cm layer of damp sphagnum moss in which the animals could construct their runways. Each arena 
contained a tubelick water bottle at 2 opposite corners; the 2 additional corners contained a food dish 
provided with the diet described above. All females constructed nests using the moss and small pices of 
cotton wadding that had been placed in one corner of each arena. Females were monitored daily at 8-hr 
intervals until parturition. 

Body mass for each neonate was recorded using a Mettler electronic balance. After birth, 
interactions between mothers and their young (in litters that had not been manipulated in any way) were 
observed daily at 2-hr intervals for a period of 22 days, and the various components of maternal behavior 
were identified. Litter sizes ranged from 3—6 (mean: 4.1). All observations were recorded manually 
through a one-way mirror to minimize disturbance to animals. 

A series of experiments were then conducted on another group of pregnant females in order to assess 
the effects of LGC on maternal behavior. In B. carolinensis, weaning occurs by 23-24 days of age 
(Punzo, unpubl. data), and the sex of the offspring can be ascertained by day 13 after birth, when the 
mammae first appear in females (personal observ.). On day 13, 13 litters with 5 pups each were culled to 4 
and randomly assigned to one of 3 treatment groups: MM (4 male young), MF (2 males and 2 females), 
and FF (4 females). Behavioral interactions between mothers and their young were observed during the 
first 3 hr of the dark phase, starting on day 2 after birth and continued every fifth day, until day 22, a day 
before the surviving young began to leave their nest. All observations were recorded manually through 
a one-way mirror to minimize disturbance to animals. The frequency of the various components of 
maternal behavior were recorded as the number and duration of episodes for each observation period. 

All statistical procedures followed those described by Sokal and Rohlf (1995). Litters were assigned 
to treatment groups randomly, a Bartlett’s test on the data showed homogeneity of variances, and a 
G-test indicated that error variances were distributed normally. As a result, 2-way analyses of variance 
(ANOVAs) considering the 3 treatment groups (LGC: MM. MF, and FF) and 7 repeated measures (age of 
litter) were used for duration and frequency of each of the behavioral components of maternal behavior. 
One-way ANOVAs (3 LGCs) were then conducted for each age. Post-hoc comparisons were performed 
using Tukey’s test. Differences between the body mass of individual neonates as well as total litter mass 
were assessed using Student t-tests. 
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REsuLTS—The body mass of neonates at birth ranged from 0.67—0.74 g (mean: 
0.69 + 0.11 SE). The mean mass for all 10 litters on day 1 was 2.82 + 0.41 g (range: 
2.43—3.17 g). There were no significant differences in mass between individual 
neonates and litters (P > 0.50). At 2 days of age, pups weighed between 1.18—1.26 g 
(mean: 1.22 + 0.07 SE). These values increased to 2.92—3.53 g (mean: 3.23 + 
0.11 g) by day 4. At 13 days of age (when pups could be sexed), pups weighed 
between 6.93—7.34 g (mean: 7.13 + 0.12), with no significant gender differences 
(P > 0.60). Weaning occurred 22—23 days after birth, and the pups weighed 8.84— 
9.47 g (mean: 9.11 + 0.23) at that time, with no significant differences between 
sexes (P > 0.50). ' 

The following behavioral components of maternal behavior were identified 
from observations on dams whose litters had not been manipulated in any way: 
nursing (young attached to nipples), licking (female licks any part of the young’s 
body), nest maintenance (any manipulation of nesting materials while young are in 
the nest), retrieval of young (female carrying any pup who had moved out of the 
central nest area back to the nest), and miscellaneous activities (these include eating, 
drinking, grooming, moving about the arena without carrying young, scratching, and 
sleeping). The pups frequently emitted high-pitched sounds, most frequently during 
the dark phase, and spent between 4—6 hr/day suckling. They slept between 14-16 
hr/day. It was during these sleeping periods that the maternal parents would leave 
their nests to feed and drink. 

During nursing bouts, the mother’s body was slightly arched over at least 1 pup, 
and the pups were attached to her nipples. Mothers directed most of their licking to 
the anogenital regions of their pups. Nest maintenance involved either carrying 
material toward the nest or manipulating materials already present in the nest. Any 
pup that moved away from the nest was grasped at the back of its neck in the 
mother’s mouth (retrieval of young) and carried back to the nest. 

Considering all observed categories of behavior for all treatment groups, 
maternal behavior was more pronounced in the days immediately following birth 
than in subsequent days (Table 1). There was a significant decrease in the frequency 
(number of episodes) of nursing (F = 24.2, P < 0.01) and licking (F= 15.3, P < 0.01) 
by maternal parents over time, and a significant increase in miscellaneous activities 
(F = 32.3, P < 0.01). Duration values also decreased significantly over time for 
nursing (F= 28.4, P < 0.01), licking (F= 12.2, P < 0.01), and nest maintenance (F = 
30.7, P < 0.01), and increased for miscellaneous activities (F = 14.2, P < 0.02). None 
of the treatment groups exhibited any retrieval of young after day 12. 

With respect to nursing, post hoc tests indicated that MM mothers exhibited 
a significantly higher frequency and duration of nursing than those from the MF and 
FF groups (P < 0.02). Frequency of licking was also greater in MM mothers that 
those from other treatment groups (P < 0.01). 

There was no significant effect of LGC on the nest maintenance activity 
or retrieval of young by maternal parents (P > 0.50). With respect to miscella- 
neous activities, post hoc tests showed that MM mothers were significantly 
more involved in non-maternal behaviors as compared to MF and FF mothers 
(P1<,.0.05), 
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TaBLe 1. Number of episodes and duration of each episode (sec) for behavioral components of 
maternal behavior observed in mothers of Blarina carolinensis rearing litters culled at birth and randomly 
assigned to one of 3 treatment (litter gender composition) groups: males only (MM—4 males); males and 
females (MF—2 of each sex); and females only (FF—4 females). Data expressed as means + SE (for 3-hr 
observation bouts over a 22-day period) for 13 litters in each treatment group. N = nursing; L = licking; 
NM = nest maintenance; RY = retrieval of young; MA = miscellaneous activities). See text for details. 


Treatment group N IL, NM RY MA 
MM 
Episodes 
Day 2 SYiS)ee Vil 62.2 + 6.9 16+ 0.4 Yai az 057) 21.6 + 6.7 
Day 7 34.4 + 7.9 45.8 + 4.6 12+ 04 17+ 0.4 BP), 22 DS) 
Day 12 24.6 + 3.7 ea 2S Aol 1.5 + 0.5 OKO 25 (01 Sys) 2S yp) 
Day 17 43" 2 16.9 + 0.8 13) a=e(08) 0 56.8 + 6.6 
Day 22 PAG =a (2 pd oe lal 1.4 + 0.3 0 10s) 22 
Duration 
Day 2 597/6-= 245 STAr = Ts /] 67.3 = 8.8 29.2 = 3.4 144.3 + 17.2 
Day 7 351.8 + 14.8 47.8 + 6.8 27.6 + 4.5 33°3° 4°45 15327 2202 
Day 12 DSS eee PASS 59338 8.4 + 0.7 MS\L8) 2S. Sol) 194.5 + 18.7 
Day 17 74S =. 21.6 DileAwets Di) Du 280),8) 0 237.8 + 26.3 
Day 22 44.5 + 8.9 Wes) 28.025) US) = ill 0 224.9 + 20.4 
MF 
Episodes 
Day 2 DN Sy Sey 39.5) = 4.9 Ove== 02 14+ 0.3 10.6 + 2.6 
Day 7 ID Grt62 Asxa]l 22. So Os) 2280), gcse O)es) 22.4 + 4.3 
Day 12 1hs3| 23.4 14.4 + 2.5 0 ele=e Ow M38) BE, 3y8) 
Day 17 koa Nes) OLD E23 0 0 34.2 2 5.1 
Day 22 2.9 + 0.4 48+ 1.8 0 0 37-6525 (6:5 
Duration 
Day 2 Sip), BE BBS) 40.2 + 5.4 DO = 7a] BS at eas 136.7 + 14.9 
Day 7 2803 205 S))8) 2m (67 34.5 + 5.5 30.4 + 7.7 148.4 + 16.8 
Day 12 190.4 + 15.8 1525922374 0 29.4 + 5.8 205.7 + 24.8 
Day 17 44.3 + 8.6 IWB8) 2 LD) 0 0 DS py 2 UY 
Day 22 BAPE BEES 14 2.3 0 0 MBE) BE DTS) 
FF 
Episodes 
Day 2 fs ape) 5308 = 8a Oloz=t0M Os OZ WESeze ES 
Day 7 21.4 + 4.7 25.1 + 4.4 0.4 + 0.1 2.3 + 0.4 DANS eee 
Day 12 | he, 2S) Test, Ecos 0 ps)| es (9 Phe\45) EAD! 
Day 17 45+04 74+ 0.9 O42 0A 0 B27] ==" 83 
Day 22 Se) 205 S28 =O) 0 0 40.3 + 6.9 
Duration 
Day 2 502.8 + 26.4 38.092 1622 SL) 22 Ooi) 28.7 + 4.8 1383/7 =s10% 
Day 7 2550 232.4 Dee) eB AY) 45.6 + 7.8 Pell 22 Bis) 155.3) 22 20,5 
Day 12 177.4 + 16.8 14.6 = 2.9 0 SOS ef 188.6 + 25.6 
Day 17 SS 853 boy 22 Bho So) li-28 Dee 0 240.3 + 32.5 
Day 22 5Gi8: 22 7-5 6.7 + 0.4 0 0 244.2 + 26.3 


DiscussioN—The results of these experiments show that maternal behavior in 
B. carolinensis depends both on the age of the pups and on LGC. Mothers of this 
shrew spend more time with younger pups and tend to male pups more frequently and 
for longer periods of time than female pups. This is the first demonstration of the 
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effect of LGC on maternal behavior for any shrew, although similar effects have been 
reported for murid rodents (Mendl and Paul, 1991; Laviola and Alleva, 1995). 
Specifically, MM litters received a greater degree of nursing and licking than MF and 
FF litters. Increased frequency of nursing may reflect the importance to overall 
survivorship of rapid growrth rates during the initial days of life, and licking is known 
to facilitate parent-offspring bonding via olfactory imprinting in numerous mammals 
(Hall, 1991). However, why female shrews lick pups from all-male litters more 
frequently than those from mixed or all-female litters is not immediately apparent. 
The bias in maternal behavior directed at all-male litters does not appear to be due to 
differences in pup development because there were no significant differences in pup 
body mass between any of the treatment groups. Thus, LGC influences maternal 
behavior in B. carolinensis but not the physical development of offspring. 

In rodents, where most of the research has been directed, it has been shown 
that male rats reared in litters containing only males exhibited impaired organization 
of adult reproductive behavior (Birk and Sadler, 1987; Hard and Larsson, 1969). In 
addition, in response to handling they defecated more frequently and emitted more 
sounds than rats reared with female littermates. Female rats reared in same-sex litters 
exhibited impairments in subsequent maternal behavior (Sharpe, 1975), while 
female mouse pups entered puberty earlier than females reared with several male 
siblings (Drickamer, 1976). 

Previous studies have shown that early social environment can have a profound 
affect on subsequent social behaviors in many mammals including shrews (Punzo 
and Chavez, 2003), mice (Laviola and Alleva, 1995), rats (Moore and Power, 1992), 
gerbils (Clark et al., 1993), cats (Caro, 1981), and primates (Lott, 1991). In addition 
to maternal behavior, these social behaviors include exploratory behavior, grooming, 
mating, aggression, play-fighting, as well as emotionality and neophobia. Future 
studies should examine the relationship between LGC and these types of behavior 
in B. carolinensis as well as other species of shrews in order to assess the degrees 
of similarity or difference between insectivores and other mammalian taxa. 


ACKNOWLEDGMENTS—I wish to thank D. Martin, Barbara Martin, and anonymous reviewers for 
comments on an earlier draft of the manuscript, L. Ludwig, C. Farmer and J. Curtis for assistance in 
maintaining the animals in captivity and conducting observations of maternal behavior, and B. Garman for 
consultation on statistical analyses. All protocols followed those described in the Gudelines for the 
Ethical Treatment of Animals (Animal Behavior Society, Washington, D.C.). 


LITERATURE CITED 


ALLEVA, E., A. CAPRIOLI, AND G. LAvIOLA. 1989. Litter gender composition affects maternal behavior of 
primiparous mouse dam (Mus musculus). J. Comp. Psychol. 103:83-—87. 

BarreTT, L., R. DUNBAR, AND P. DuNBAR. 1992. Environmental influences on play behavior in immature 
gelada baboons. Anim. Behav. 44:111—115. 

Birk, L. AND D. SADLER. 1987. Differences in maternal behavior of rats and the sociosexual development 
of the offspring. Develop. Psychobiol. 20:85—99. 

Caro, T. M. 1981. Sex differences in the termination of social play in cats. Anim. Behav. 29:271-—279. 

CHURCHFIELD, S. 1990. The Natural History of Shrews. Cornell Univ. Press, Ithaca, NY, ,178 pp. 

CuarK, M. M., P. KarPLuK, AND B. G. GALEF. 1993. Hormonally mediated inheritance of acquired 
characteristics in Mongolian gerbils. Nature 364:712—713. 


No. 4 2003] PUNZO—SHREW BEHAVIOR Sls 


De Caso, C., M. Ketty, M. KITCHEN, AND P. Viveros. 1992. Effects of neonatal sexual isolation on 
u-receptor development and nociceptive responses in the rat. Neurosci. Res. Comm. 10:79-86. 

DENENBERG, V. H. AND J. R. Morton. 1964. Infantile stimulation, prepubertal sexual-social interactions, 
and emotionality. Anim. Behav. 12:11-13. 

DRICKAMER, L. C. 1976. Effect of size and sex ratio of litter on the sexual maturation of female mice. 
J. Reprod. Fertil. 46:369-374. 

Genoways, H. H. Anp J. R. CHoaTe. 1998. Natural history of the southern short-tailed shrew, Blarina 
carolinensis. Occ. Pap. Mus. Southwest. Biol. 8:1-43. 

GiNGERICH, J. L. 1999. Florida’s Fabulous Mammals. World Publs., Tampa, FL, 128 pp. 

HALL, G. 1991. Perceptual and Associative Learning. Oxford Univ. Press, Oxford, 412 pp. 

Harp, E. AND K. Larsson. 1969. Dependence of adult mating behavior in male rats on the presence of 
littermates in infancy. Brain Behav. Evol. 1:405-419. 

LavIoLA, G. AND E. ALLEvA. 1995. Sibling effects on the behavior of infant mouse litters. J. Comp. 
Psychol. 109:68—75. 

Lott, D. F. 1991. Intraspecific Variation in the Social Systems of Wild Vertebrates. Cambridge Univ. 
Press, New York, NY, 388 pp. 

MENDL, M. AND E. S. PauL. 1991. Litter composition affects parental care, offspring growth, and the 
development of aggressive behavior in wild house mice. Behaviour 116:90—109. 

Moore, C. L. AND K. L. Power. 1992. Variation in maternal care and individual differences in play, 
exploration, and grooming of juvenile Norway rat offspring. Develop. Psychobiol. 25:165—182. 

O’FarrRELL, M. J., D. W. KAUFMAN, J. B. GENTRY, AND M. H. Smitu. 1977. Reproductive patterns of some 
small mammals from South Carolina. Florida Scient. 40:76-84. 

Punzo, F. 2003. Behavioral interactions between the shrew Blarina carolinensis and four species of 
Florida mice in the presence and absence of mouse nestlings. Florida Scient. 66:63—67. 

PuNZzo, F. AND S. CHAvez. 2003. Effect of aging on a spatial learning task and running speed in the least 
shrew, Cryptotis parva. J. Mammal. (in press). 

SHARPE, R. M. 1975. The influence of the sex of littr mates on subsequent maternal behaviour in Rattus 
norvegicus. Anim. Behav. 23:551—-559. 

SOKAL, R. R. AND F. J. RouiF. 1995. Biometrics. 3rd. ed. W. H. Freeman, San Francisco, CA, 818 pp. 

WhirTAKER, J. O., JR. 1996. Field Guide to North American Mammals. Alfred A. Knopf, NY, 938 pp. 


Florida Scient. 66(4): 308-313. 2003 
Accepted: March 28, 2003 


Environmental and Chemical Sciences 


SEDIMENT CONTAMINATION SURVEY ON ST. MARKS 
NATIONAL WILDLIFE REFUGE 


JON M. HEMMING* AND MICHAEL S. BRIM 


United States Fish and Wildlife Service, Panama City Field Office, 1601 Balboa Avenue, 
Panama City, FL 32405 


ABSTRACT: A survey was conducted by the U.S. Fish and Wildlife Service to assess habitat quality 
on the St. Marks National Wildlife Refuge (NWR). Sediment samples were collected at 32 sites within 
(n= 14) or north (n= 18) of the NWR. The survey showed only site and use-specific contamination. Much 
of the surveyed area possessed little or no contaminant residues. Sediment composition around St. Marks 
NWR was comparable to coastal areas of the northeastern Gulf of Mexico. Metal contamination of 
sediments was not found on-refuge but was found at 6 off-refuge sites and included moderate 
concentrations of copper and mercury. Polycyclic aromatic hydrocarbon (PAH) contamination was also 
not found on-refuge, however, off-refuge PAH contamination was slightly more widespread than metal 
contamination being found at 9 off-refuge sites. No organochlorine contamination was detected in 
samples taken on or off-refuge. Sixteen sites, including both on and off-refuge sites, were found to have 
relatively high concentrations of aliphatic hydrocarbons. Oil and grease contamination was found at 11 
sites total, but at only I site on the NWR. The survey objective was to provide baseline information from 
which to determine the need for additional monitoring and for use in developing management strategies. 


Key Words: St. Marks National Wildlife Refuge, sediment, contamination, copper, 
mercury, PAH 


THE St. Marks National Wildlife Refuge (NWR) was established in 1931 to 
provide wintering grounds for migratory bird species. The refuge encompasses 
nearly 68,000 acres of upland habitat in Jefferson, Wakulla and Taylor Counties, as 
well as 32,000 aquatic habitat acres in Apalachee Bay (Figure 1). St. Marks NWR is 
one of the oldest in the National Wildlife Refuge System. In addition to providing 
habitat for migratory birds, St. Marks is home to diverse plant and animal 
communities and Federal and State listed Threatened and Endangered plants and 
animals. St. Marks NWR inhabitants include southern bald eagle (Haliaeetus 
leucocephalus), red-cockaded woodpecker (Picoides borealis), least tern (Sterna 
antillarum), woodstork (Mycteria americana), American alligator (Alligator 
mississippiensis), eastern indigo snake (Drymarchon corais couperi), swallow- 
tailed kite (Elanoides forficatus), peregrine falcon (Falco peregrinus tundrius)and 
Florida black bear (Ursus americanus floridanus). The overall refuge goals include: 
providing winter habitat for migratory birds and waterfowl, habitat for endangered 
species, and habitat for all of its resident wildlife. 

Managing habitat quality on St. Marks NWR is essential to the mission of the 
U.S. Fish and Wildlife Service (Service) and the goals of the refuge. Through 
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Fic. 1. St. Marks National Wildlife Refuge (outlined) on nearly 68,000 upland acres in Jefferson, 
Wakulla, and Taylor Counties, as well as 32,000 acres of Apalachee Bay. 


management designed to promote ecological integrity, the Service strives to protect 
the fish, wildlife and habitat entrusted to it. For this reason, habitat quality surveys 
are conducted to report the status of these systems. This report describes a general 
sediment quality survey conducted in 1988 to reveal environmental contamination 
of the aquatic systems of St. Marks NWR. To this end, the Service collected 32 
sediment samples on and north of St. Marks NWR and had them analyzed for metal, 
hydrocarbon and organochlorine environmental contaminants. This survey was 
intended to elucidate the contamination status of the St. Marks NWR for the pur- 
pose of providing managers information needed to contend with issues challenging 
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TABLE 1. Sample information for sediment samples taken by the U.S. Fish and Wildlife Service on 
St. Marks National Wildlife Refuge in 1988: sample identification, sampling location description, latitude 


and longitude (degrees, minutes, hundredths of minute). 


Sample ID Location Latitude Longitude 
SM1 East River Pool 30°12.20 84°08.50 
SM2 Mounds Pond #1 30°10.17 84°09.10 
SM3 Tower Pond 30°08.83 84°09.01 
SM4 Mounds Pond #3 30°09.75 84°08.30 
SM5 East River 30°10.87 84°09.82 
SM6 East River below Denham 30°10.48 84°10.10 
SM7 Boat Basin at lighthouse 30°08.42 84°10.42 
SM8 Stony Bayou #2 30°11.33 84°07.10 
SM9 Stony Bayou #1 30°12.17 84°08.30 
SM10 Picnic Pond 30°05.22 84°09.80 
SM11 St Marks River at Buoy 27E 30°06.25 84°11.32 
SM12 St Marks River at Buoy 27W 30°06.18 84°11.55 
SM13 St Marks River at Buoy 42E 30°07.63 84°11.81 
SM14 St Marks River at Oliver Bayou 30°07.63 84°11.81 
SM15 Big Boggy Bayou 30°10.30 84°13.06 
SM16 Wakulla River below SMYC* 30°09.59 84°12.38 
SM17 Wakulla River at Shell Island 30°09.83 84°12.67 
SM18 Wakulla River 2 miles from mouth 30°10.43 84°13.74 
SM19 St Marks River above power plant 30°10.30 84°11.09 
SM20 St Marks River at East Side 30°09.49 84°12.10 
SM21 St Marks River at fuel docks 30°09.60 84°11.97 
SM22 St Marks River at Canal Marina 30°09.74 84°11.69 
SM23 St Marks River at turning basin S 30°09.68 84°11.61 
SM24 St Marks River at marina 30°09.73 84°11.60 
SM25 St Marks River at turning basin N 30°09.72 84°11.57 
SM26 St Marks River at new marina N 30°09.82 84°11.63 
SM27 St Marks River at fuel loading 30°09.95 84°11.58 
SM28 Power Plant at discharge 30°10.02 84°11.60 
SM29 Power Plant at fuel load 30°10.10 84°11.53 
SM30 Power Plant at intake 30°10.25 84°11.58 
SM31 St Marks River below Newport Bridge 30°12.19 84°10.38 
SM32 St Marks above Newport Bridge 3012253 84°10.25 


* St. Marks Yacht Club. 


habitat quality. Secondarily, the gathered data were intended to be used as baseline 
for future evaluations and assessments of St. Marks NWR system health. 

Habitat quality of aquatic systems has been evaluated via surveys of sedi- 
ment contamination (O’Connor, 1991; USNOAA, 1991; Bolton et al., 1985). The 
challenge lies in the interpretation of ecological risk posed by sediment contami- 
nation (Long et al., 1995). Many have provided numeric guidelines based on re- 
ported effects of exposure as a means to estimate relative risk to living organisms 
(Buchman, 1999; Long et al., 1995; MacDonald, 1993; Persaud et al., 1992; Di Toro 
et al., 1991; Long and Morgan, 1990; US EPA, 1989). 

The results of this survey were compared to the findings of Long and co- 
workers (1995) to estimate risk to living resources from exposure to contaminated 
sediments and to assess overall aquatic ecosystem health on the St. Marks NWR. 
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TABLE 2. Chemical analytes measured in sediment samples taken on St. Marks National Wildlife 


Refuge 1988. 
Polycyclic Aromatic Aliphatic Organochlorines 

Metals hydrocarbons hydrocarbons and pesticides 
*Silver *Napthalene n-Dodecane Hexachlorobenzene 
Aluminum *Fluorene n-Tridecane a,b,g and d-BHC 
* Arsenic *Phenanathrene n-Tetradecane Oxychlordane 
Boron * Anthracene Cyclohexane Heptachlor 
Barium *Fluoranthrene Pentadecane a, g-Chlordane 
Beryllium *Pyrene Cyclohexane t-Nonachlor 
*Cadmium *Benz(a)anthracene n-Hexadecane Toxaphene 
*Chromium *Chrysene n-Heptadecane *Total PCBs 
*Copper Benzo(b)fluoranthrene Pristane *DDT analytes 
Iron Benzo(k)fluoranthrene n-Octadecane Dieldrin 
*Mercury Benzo(e)pyrene Phytane Endrin 
Magnesium *Benzo(a)pyrene n-Nonadecane cis-Nonachlor 
Manganese Dibenzo(a,h)anthracene N ecosane Mirex 
Molybdenum Benzo(g,h,i)perylene Total Ahs Dicofol 
*Nickel *Total PAHs Dicamba 
*Lead Dichloprop 
Selenium Silvex 
Strontium 2,4-D 
Thallium 2,4,5-T 
Vanadium 2,4-DB 
*Zinc Pentachlorophenol 


* Sediment Quality Guidelines available from Long et al. 1995. 


Long and co-workers (1995) developed Effects Range Low (ERL) and Effects 
Range Median (ERM) guidelines for evaluating sediment contamination. Sediment 
contaminant concentrations exceeding the criterion ERL indicated that adverse 
negative effects on living resources may increase in incidence from rare to 
occasional. Sediment contaminant concentrations exceeding the ERM may indicate 
that adverse effects will likely occur. 


MeETHops—Sediment sampling was chosen to evaluate habitat quality on St. Marks NWR because 
of the tendency of many environmental contaminants to accumulate in sediments, thus providing an 
established route of entry into system food webs. Thirty-two sediment samples were collected at various 
sites from waters on (n= 14) and north (n= 18) of St. Marks NWR (Table 1). Sites 1- to 14 were collected 
on St. Mark NWR and sites 15—32 were collected in the more industrialized area of St. Marks River north of 
the NWR. Sediment samples were composite 1-liter samples consisting of five 200 ml subsamples. 
Samples were collected using a standard ponar 316 stainless steel grab for all samples except those in water 
too shallow for the contaminant survey boat. Samples not taken with the standard ponar were taken with 
a petite ponar stainless steel grab. Depth of sediment samples collected depended on the ponar used and 
type of sediment at each station, but all were less than 10 cm. Samples collected in the field were 
immediately put into laboratory-certified, chemically-cleaned, 1-liter amber, glass jars with Teflon-lined 
lids, and placed on ice in coolers. Samples were temporarily stored at the Panama City Field Office in 
freezers at —23° C until shipment to analytical laboratories. Sediment samples were analyzed for metal, 
hydrocarbon (aliphatic and aromatic), and organochlorine chemical contaminants by the Geochemical and 
Environmental Research Group at Texas A&M University. A complete list of analytes is provided in Table 
2. Additional samples were analyzed for particle size and total organic carbon (TOC). 
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TABLE 3. Specific information on sites with metal analytes equal to or exceeding recommended 
sediment quality guidelines (Long et al., 1995): site, analyte, sediment concentration (ug/kg, ppb, dry 
weight), effects range low (ERL), and effects range median (ERM). 


Site Analyte Sediment Concentration ERL ERM 
SM28 Copper 39.3 34.0 270 
SM24 Mercury 0.16 0.15 0.71 
SM25 - 0.16 3 9 
SM27 ss 0.15 ” 
SM30 3 0.16 » 


SM31 5 0.17 » » 


Data were compared to the Effects Range Low (ERL) and Effects Range Median (ERM) guidelines 
of Long and co-workers (1995) to estimate potential risk to living resources from exposure to 
contaminated sediments and to assess the overall ecosystem health of the St. Marks NWR. 

In the case of aliphatic hydrocarbon (AH) and oil and grease sediment contamination, sediment 
quality guidelines were not available. For the purpose of contamination evaluation, areas/sites were 
recommended for further evaluation based on relative sediment concentrations. If concentrations were 
moderately higher than other sites (1,000 + ppb) or higher (2,000 + ppb) than other sites surveyed, further 
evaluation was recommended accordingly. 


RESULTS AND DiscussioN—Tables are provided to show specific analytes that 
may constitute a problem in given areas. The areas of possible concern are further 
described by site, analyte, concentration, effects range low (ERL), and effects range 
median (ERM). Sediment concentrations equal to or exceeding the guidelines pro- 
vided in the Long and co-workers (1995) publication indicate occasional (ERL) or fre- 
quent (ERM) association with adverse effects on living resources. Tables containing 
ERL and ERM guidelines are presented for metals (Table 3) and polycyclic aromatic 
hydrocarbons (Table 4), but not for organochlorines because none were detected. 

Areas with relatively high concentrations of aliphatic hydrocarbons or oil and 
grease were divided into groups possibly needing further investigation, but not 
necessarily indicating risk from exposure. Risk assessment for these contaminant 
categories was omitted because there are no sediment guidelines currently available 
for these classes of chemicals. Areas recommended for further evaluation were 
based on sediment concentrations that were high relative to other sites in the survey. 
Tables summarizing these sites are shown for sediment aliphatic hydrocarbons 
(Table 5) and oil and grease (Table 6). 

Sediment composition on and around St. Marks NWR was comparable to other 
coastal areas of the northeastern Gulf of Mexico (Hemming et al., 2002; Brim et al., 
2000; Brim, 1998). The typical sandy sediments of the Gulf coast were evident with 
sand fractions measuring as high as 95%. Sand fractions were also found to be as low 
as 15%, but most often above 65%. The silt and clay fractions were similarly variable 
and reciprocal to the sand fraction. Gravel was not uncommon in samples, but 
constituted a maximum of 5.7% of the total sediment composition by weight. Per- 
cent total organic carbon ranged from 0.2 to 4.9 which was also typical of sediments 
from the northern Gulf of Mexico (Hemming et al., 2002; Brim et al., 2000; Brim, 
1998). 
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TaBLE 4. Specific information on sites with polycyclic aromatic hydrocarbon analytes equal to 
or exceeding recommended sediment quality guidelines (Long et al., 1995): site, analyte, sediment 
concentration (i1g/kg, ppb, dry weight), effect range low (ERL; pg/kg, ppb, dry weight), and effects range 
median (ERM; g/kg, ppb, dry weight). 


Site Analyte Sediment Concentration ERL ERM 
SM21 Phenanthrene 292 240 1,500 
ra Anthracene 643 85 1,100 

Ea benz(a)anthracene 468 261 1,600 

w chrysene Sia 384 2,800 

e Total PAHS’ 4,269 1,700 9,600 
SM22 Total PAHS! 1,992 1,700 9,600 
SM23 Total PAHS! 3,583 1,700 9,600 
SM24 Total PAHS! 2,190 1,700 9,600 
SM25 Total PAHS! 2,500 1,700 9,600 
SM26 Total PAHS’ 2,565 1,700 9,600 
SM27 Fluorene 429 19 540 
” Phenanthrene 1,786 240 1,500 

Me Anthracene 35)7/ 85 1,100 

“3 Fluoranthrene 1,643 600 5,100 

.. Pyrene 1,286 665 2,600 
Chrysene 500 384 2,800 

as Total PAHs' 8,571 1,700 9,600 
SM28 Total PAHs' 2,208 1,700 9,600 
SM29 Chrysene 399 384 2,800 
Total PAHs' 2,822 1,700 9,600 


" Total sum of analyzed PAHs. 


Metal contamination was found in sediment samples taken just north of St. 
Marks NWR, but not on the actual refuge. Only 2 metals, copper and mercury, were 
found at concentrations expected to increase the incidence of adverse negative 
effects on living resources from rare to occasional (Effects Range Low, ERL, Long 
et al., 1995). Copper exceeded the ERL at 1 site and mercury exceeded the ERL at 
5 sites. No metals exceeded the Effects Range Median (ERM), therefore frequent 
adverse effects were not expected from metal exposure. Metal concentrations in 
samples from all 6 sites were only slightly higher than the respective ERLs. 

Polycyclic aromatic hydrocarbon (PAH) contamination was not found on 
the refuge, but was somewhat more widespread than metal contamination at sites 
north of the refuge (9 sites). Many of the PAH contaminated sites were the same as 
those with slight metal contamination. The PAH analytes chrysene, phenanthrene, 
and anthracene were most frequently in excess of the ERL, however, they were only 
found 3, 2, and 2 times each, respectively. Phenanthrene was detected at a concen- 
tration in excess of the ERM at a fuel loading area. Contamination in excess of the 
ERM may increase the incidence of adverse effects to frequent (Long et al., 1995). 
The most common sediment guideline exceeded was the ERL guidance value for sum 
total PAHs. This criterion was exceeded at all 9 sites where PAH contamination was 
found, even at sites where no individual analytes exceeded the sediment quality 
guidelines. PAH contamination at these sites was not surprising in the industrial area 
on St. Marks River and with fuel docks and marinas nearby. PAHs are fuel fractions 
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TABLE 5. Specific information on sites with relatively high concentrations of total aliphatic 
hydrocarbon analytes for the St. Marks National Wildlife Refuge sediment quality survey, 1988: site, 
analyte, sediment concentration (g/kg, ppb, dry weight), possible need for further evaluation (possible; 
ug/kg, ppb, dry weight),and probable need for further evaluation (probable; ug/kg, ppb, dry weight). 


Site Analyte Sediment Concentration Possible Probable 
SM2 Sum Total Aliphatic 3,770 1,000 2,000 
SM3 Hydrocarbons 1,166 on ae 
SM4 . 1,306 a a 
SM7 = 3,842 55 » 
SM8 ce 4,056 6 ” 
SM12 5 1,073 59 
SM21 a 2,423 es _ 
SM22 _ 1,504 5 » 
SM23 * 1,398 5 » 
SM24 a 2,470 
SM25 es 29 x » 
SM26 5 1,832 55 » 
SM27 as DeSO\0) 55 » 
SM29 - 1,149 55 » 
SM30 a 2,475 — s 


SM31 2 1,141 » » 


and would be expected to be found where fossil fuels are burned and may enter the 
environment. 

Not only were organochlorine (OC) compounds not found to exceed the sedi- 
ment quality guidelines, but no OC contamination was detected in samples taken. 
The significance of the absence of OC residues in the sediments on and around 
St. Marks NWR is unclear. The reported historical use of OCs in pest management 
in the area makes the lack of detectable analytes noteworthy. 

AHs and oil and grease residues were the only contaminant groups found in 
sediment samples taken from the aquatic locations actually on St. Marks NWR. 
Additional contamination was found at sites north of the refuge where metal and 
PAH contamination was found. Sixteen sites were recommended for further inves- 
tigation due to the AH concentration of the sediments. Of those 16 sites, sediment 
AH contamination of 9 of the sites was considered moderately high, versus high at 
the remaining 7 sites. Like PAHs, AHs are fuel fractions and often occur where fuel 
is used or dispensed. More specifically, AHs are the lighter fuel fractions typically 
dominating small engine fuels. The co-occurrence of both PAHs and AHs at sites 
around fueling dock and marinas may be expected. AH contamination of other sites, 
including bayous, lakes and ponds, may have been the result of motorboat traffic due 
to the inefficient use of gasoline products by these 2-cycle engines. Use of the more 
recently available 4-cycle engines may help to lessen the level of contamination at 
these sites/areas. 

Oil and grease contamination, based on the above described benchmarks, was 
for the most part found around boat activity and fueling stations as described for 
PAHs. However, this type of contamination was found at a coastal bayou site on 
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TABLE 6. Specific information on sites with relatively high concentrations of combined oil and 
grease analytes for the St. Marks National Wildlife Refuge sediment quality survey 1988: site, analyte, 
sediment concentration (t1g/kg, ppb, dry weight), possible need for further evaluation (possible; pg/kg, 
ppb, dry weight), and probable need for further evalutation (probable; pg/kg, ppb, dry weight). 


Site Analyte Sediment Concentration Possible Probable 
SM8 Combined Oil 1,469 1,000 2,000 
SM18 and Grease 1,738 sy Fe 
SM21 3 1,289 i. 3 
SM22 - 1,278 i 55 
SM23 a 1,505 *) 
SM24 s 2,530 99 
SM25 % 1,880 - 99 
SM26 5 E22 5 - 
SM27 % 95 . > 
SM29 ss 1,522 5 99 


SM30 8 Pa NI 22 » » 


St. Marks NWR. Oil and grease use in lubrication of motorboat engines likely 
contributed largely to this specific contamination. 


CoNncLusIons—The sediment quality survey of the aquatic habitats of St. Marks 
National Wildlife Refuge (NWR) documented contamination of specific sites or area 
only. Much of the area surveyed possessed little or no contamination. Aliphatic 
hydrocarbons and oil and grease residues were the only contamination found in 
sediment samples taken from the aquatic environments on St. Marks NWR. 
However, contamination was found at sites in the more industrialized area of St. 
Marks River north of the NWR. 

These data provide a benchmark for subsequent investigations and activities. A 
re-evaluation of all sites on St. Marks NWR for recent changes and temporal trends 
in sediment contamination and habitat quality should be performed. Of particular 
interest, the biological impacts of identified contaminated sediments should be 
evaluated. Evaluations should also assess potential adverse ecological effects from 
exposure to aliphatic hydrocarbon and oil and grease residues in sediments. These 
data can be used further in the development of sediment quality guidelines for classes 
of chemicals for which guidelines do not exist. 
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